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RESEARCH  ON  NONDESTRUCTIVE  TESTING 


i 


I.  INTRODUCTION 

The  aim  of  this  year's  work  has  been  to  establish  new  techniques  for 
measuring  the  properties  of  materials  with  the  ultimate  aim  of  determining 
material  properties  such  as  hardness  and  yield  strength  and  being  able  to 
predict  the  life  and  servicibil ity  of  mechanical  structures.  We  have 
limited  ourselves  to  acoustic  measurement  techniques.  Our  purpose  is  to 
correlate  the  results  of  acoustic  measurements  with  the  mechanical  proper¬ 
ties  of  the  material  of  interest. 

The  project  has  been  divided  into  two  parts:  (1)  the  characterization 
of  microstructure  by  acoustic  NDE,  and  (2)  the  measurement  of  stress 


states  by  acoustic  NOE. 


The  first  year's  work! 
very  successful.  We  have 


on  characterization  of  microstructure  has  been 
shown  that: 


1.  Relative  velocity  change  measurements  can  be  used  to  map  micro- 
structure  with  precision  in  a  given  piece  of  steel.  We  have  made 
the  first  measurements  of  hardness  variation  along  a  Jominy  test 
bar  by  acoustic  methods.  These  results  are  in  excellent  agreement 
with  the  more  conventional,  much  slower  technique  using  a  Knoop 
hardness  indentor  for  determining  hardness. 

2.  We  have  shown  that  absolute  velocity  measurements  can  be  carried 

out  extremely  accurately.  But  they  are  normally  subject  to  large 

random  variations  from  sample  to  sample,  due  to  minor  variations 

of  the  constituents  of  the  alloys  in  typical  structural  steels. 
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3.  Acoustic  attenuation  measurements  are  very  sensitive  to  micro¬ 
structure  and  have  good  potential  for  practical  microstructural 
NOE.  We  have  developed  an  automatic  method  for  measuring 
attenuation  in  materials  over  a  wide  frequency  range  which 
ultimately  should  enable  us  to  measure  microstructural  variations 
through  the  depth  of  a  metallic  sample.  This  should  have  appli¬ 
cations  to  such  measurements  as  those  of  case  hardening. 

4.  During  the  course  of  development  of  these  measurements,  we  have 
improved  our  existing  computer  controlled  velocity  measurement 
system  so  that  it  can  routinely  measure  relative  velocity  to 

1  part,  in  10^  at  any  point  of  a  flat  metal  sample. 

Our  techniques  for  measuring  the  stress  state  of  the  materials  have 
also  made  good  progress.  We  have  developed  a  fairly  complete  theory 
which  shows  that  it  is  possible  to  measure  by  acoustic  methods  the 
so-called  J  and  M  energy  integrals  in  the  neighborhood  of  a  crack.  The 
evaluation  of  these  integrals  makes  it  possible  to  determine  whether  a 
crack  will  rapidly  grow  when  a  given  stress  is  applied  to  a  structure,  a 
property  which  is  normally  very  difficult  to  measure.  In  related  experi¬ 
mental  work  on  another  contract,  we  have  carried  out  experiments  which 
have  confirmed  some  of  our  initial  hypotheses  and  which  show  that  indeed 
it  is  possible  to  measure  these  J  and  M  integrals.  On  this  grant,  we 
developed  the  basic  theory  required  for  this  purpose. 

In  addition,  we  have  developed  a  completely  new  theory  which  gener¬ 
alizes  these  concepts  to  piezoelectric  materials  and  to  thermoelastic 
effects  and  magnetoelastic  effects.  Some  unexpected  results  have  been 
obtained  wnich  in  fact  turned  out  to  have  been  measured  experimentally 
some  time  ago  in  unpublished  work.  For  instance,  it  was  shown  that  the 
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presence  of  an  electric  field  in  a  piezoelectric  material  could  affect 
the  rate  of  growth  of  a  crack. 

The  theories  developed  are  far  more  general  than  are  required  for 
the  evaluation  of  the  single  crack.  As  we  have  shown,  it  is  also  pos¬ 
sible  to  determine  whether  cracks  are  present,  and  to  determine  the 
energy  associated  with  a  large  number  of  microcracks  in  a  fatigued 
structure  or  a  welded  structure. 

Thus,  in  summary,  we  are  in  the  process  of  developing  entirely  new 
measurement  techniques  for  metal,  ceramic,  and  other  solid  structures. 


II.  CHARACTERIZATION  OF  MICROSTRUCTURE  RY  ACOUSTIC  NDE 
(J.  Shyne,  6.  S.  Kino,  D.  11  ic,  F.  Stanke,  N.  Grayeli) 


A.  Introduction 

Mechanical  properties  of  materials,  such  as  tensile  strength,  hard¬ 
ness,  and  impact  toughness,  are  strongly  dependent  on  microstructural 
features  like  grain  size  and  shape,  the  proportions  and  spatial  distri¬ 
bution  of  phases  present,  and  macroscale  heterogeneities  in  microstruc¬ 
ture  resulting  from  local  compositional  or  t henna  1 -mechanical  processing 
variations.  This  is  true  of  both  metallic  and  ceramic  materials.  The 
microstructural  character  of  materials  and  their-  related  properties  are 
usually  assessed  by  microscopic  examination  or  by  direct  measurement  of 
properties;  such  procedures  are  often  costly  and  are  necessarily  destruc¬ 
tive,  requiring  sacrificial  example  specimens. 

We  have  made  considerable  progress  toward  the  goal  of  using  acoustic 
techniques  to  characterize  microstructure  in  materials  nondestructively. 
The  greatest  success  has  been  gained  by  means  of  acoustic  velocity 
measurements;  however,  acoustic  attenuation  measurements  are  also  being 
used.  In  the  long  run,  attenuation  measurements  ought  to  be  more  versatil 
and  more  useful  than  velocity  measurements.  Progress  to  date  with 
attenuation  measurements  has  mostly  been  in  equipment  and  technique 
development. 


R^  Measurement  Techniques 

Significant  improvements  and  new  developments  have  been  introduced  in 
our  techniques  for  very  precise  measurements  of  velocity  and  attenuation 


of  acoustic  waves  propagating  through  solid  samples.  All  these  measure¬ 
ments  are  performed  using  commercial  non-bonded  transducers,  using 
multiple  pulse-echo  methods,  so  that  the  effects  of  the  liquid  buffer 
are  cancelled  out  from  the  measurements.  The  measurements  are  also 
interfaced  to  an  on-line  computer,  which  performs  data  collections, 
reduction,  and  correction  for  extraneous  effects  (e.g.  diffraction, 
transducer  response,  buffer/sample  reflections),  and  presentation  of  two- 
dimensional  velocity  fields  in  the  sample,  or  frequency  dependence  of  the 
corrected  attenuation. 

1.  Velocity  Measurement  Technique 

The  electronic  signal  processing  section  of  our  existing  set-up  was 
completely  rebuilt  in  order  to  eliminate  problems  with  long-term  drifts 
which  had  limited  its  accuracy  in  the  past.  In  particular,  the  receiver 
chain  was  completely  replaced  with  a  variable-gain  amplifier  which  enables 
us  to  measure  velocity  in  samples  with  wide  variation  in  losses;  the 
detector  section  was  rebuilt  to  improve  the  rejection  of  high-frequency 
noise;  a  lock-in  amplifier  with  a  tuned  signal  channel  replaced  the  old 
combination  of  a  linear  gate  and  tuned  amplifier,  and  an  analog  phase- 
lock  loop  was  introduced  at  the  system  output.  All  of  these  changes  have 
resulted  in  more  than  an  order  of  magnitude  improvement  of  the  system 
precision,  to  the  present  value  of  10"^.  This  is  much  smaller  than  the 
effects  which  we  are  measuring  in  both  our  studies  of  stress  and  micro¬ 
structure,  which  are  of  the  order  of  0.1%  or  larger. 

2.  Absolute  Velocity  Measurements 

Five  specimens  of  plain  carbon  steel  were  prepared  by  machining  1  cm 
thick  flat  plates  oriented  both  longitudinally  and  transversely  relative 


to  the  steel  bar  stock.  Since  the  acoustic  wave  propagation  direction 
was  through  the  1  cm  plate  thickness,  acoustic  velocity  could  be  measured 
both  parallel  to  and  transverse  to  the  rolling  direction  for  each  steel 
composition  to  detect  the  influence  of  any  preferred  orientation  or 
crystal  texturing.  The  compositions  of  the  steels  are  shown  in  Table  I. 
All  specimens  were  given  identical  heat  treatment;  they  were  heated  to 
900°C  and  air  cooled  (normalized).  This  resulted  in  homogeneous  pearl ite/ 
ferrite  microstructures,  the  pearl ite  varying  from  12  to  100  volume  pet. 
over  the  range  of  carbone  contents.  After  machining,  the  flat  specimens 
were  lapped  to  assure  that  their  flat  surfaces  were  parallel  to  within 
2.5  pm. 

The  resulting  velocity  measurements,  which  include  correction  for 
diffraction  effects,  exhibited  random  scatter  that  total ly  obscured  the 
variation  of  velocity  anticipated  from  the  variation  in  microstructure. 

One  can  easily  calculate  an  expected  acoustic  velocity  from  the  known 
densities  of  pure  iron  ferrite  and  Fe^C,  and  known  elastic  constant 
data  for  ferrite/Fe-jC  mixtures.^  Increasing  carbon  should  cause  a 
decrease  in  the  longitudinal  acoustic  velocity  linearly  proportional  to 
carbon  content;  velocity  should  decrease  0.85  pet.  per  wt.  pet.  carbon 
in  the  steel,  all  other  factors  being  identical.  The  scatter  observed 
in  the  measured  velocity  was  not  caused  by  variations  in  preferred 
orientation.  Had  there  been  any  significant  degree  of  preferred  orienta¬ 
tion,  the  longitudinal  and  transverse  velocities  would  differ  much 
more  than  the  slight  variations  observed.  We  believe  that  the  random 
variations  in  such  elements  as  manganese,  silicon,  sulfur,  and  others 
always  present  in  at  least  trace  amounts,  cause  variations  in  density 
and  elastic  modulus  and  thus  introduce  random  variations  in  acoustic 
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velocity  equal  to  or  greater  than  the  systematic  effect  of  carbon  in 
changing  the  microstructure.  For  example,  we  have  calculated  that  the 
range  of  manganese  contents  in  our  specimens  (see  Table  I)  causes 
density  variations  large  enough  to  obscure  all  the  above  change  attrib¬ 
utable  to  variations  in  microstructure. 

Therefore  we  conclude  that  absolute  velocity  measurements  are  not 
practical  means  for  characterizing  steel  microstructure,  because  random 
compositional  variations  will  confuse  any  attempt  to  calibrate  absolute 
acoustic  velocity  variations  with  microstructure.  Similar  small  compo¬ 
sitional  differences  also  occur  in  nonferrous  alloys;  this  can  be  expected 
to  render  absolute  velocity  measurements  equally  impractical  for  the 
microstructural  characterization  for  such  materials  as  is  the  case  for 
steels. 


TABLE  I 

Composition  of  Plain  Carbon  Steels  Used  For 
Acoustic  Velocity  Measurements 


Steel  Type 
(AISI ) 

C 

Mn 

Si 

P 

S 

1010 

0.10 

0.45 

0.03 

0.007 

0.032 

1020 

0.20 

0.52 

0.21 

0.009 

0.028 

1035 

0.33 

0.72 

0.17 

0.018 

0.025 

1060 

0.55 

0.80 

0.19 

0.012 

0.023 

1095 

0.93 

0.50 

0.22 

0.009 

0.030 
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3.  Relative  Velocity  Measurements 

If  the  material  composition  is  uniform  within  the  object  (as  is 
usually  the  case),  variations  in  acoustic  velocity  can  be  unequivocally 
related  to  microstructural  variations  from  one  location  to  another. 

There  are  many  situations  where  microstructure  is  not  uniform  and  the 
pattern  of  microstructural  variation  is  important,  e.g.  in  an  induction 
hardened  steel  shaft  heat  treated  to  be  martensitic  on  the  outer  surface 
but  pearl itic  in  its  interior.  Jominy  end-quench  test  bars  provide 
another  example  of  steel  objects  containing  microstructural  gradients. 

We  successfully  used  relative  velocity  measurements  to  survey  quantita¬ 
tively  microstructure  of  end-quench  test  bars  of  several  different  alloy 
steel  compositions. 

The  Jominy  Test  or  End-Quench  Test  is  the  standard  metallurgical 
quality  control  procedure  for  measuring  the  hardening  response  of  heat 
treatable  steels.  The  end-quench  test  specimen  is  a  one- inch  diameter 
round  bar,  four  inches  long.  The  bar  is  heated  to  about  E  ’C.  It  is 
then  placed  in  a  fixture  and  cooled  by  a  jet  of  cold  water  impinged  upon 
one  end.  This  results  in  uniaxial  heat  flow  toward  the  water-quenched 
end,  and  reproducible  cooling  rates  that  decrease  with  increasing  distance 
away  from  the  rapidly-cooled  quenched  end.  The  fast  cooling  rate  at  the 
quenched  end  causes  the  formation  of  hard,  strong  martensite;  the  slower 
cooled  end  transforms  to  softer,  weaker  pearl ite  or  a  pearl ite  and  ferrite 
mixture,  depending  on  composition.  At  intermediate  locations  mixed 
martenslte/pearl ite  microstructures  result  from  the  intermediate  cooling 
rates.  The  position  of  the  transition  from  martensite  to  pearl ite  is  a 
measure  of  the  hardening  response.  The  usual  way  of  assessing  the  micro- 
structural  gradient  along  the  length  of  an  end-quench  bar  is  to  survey 
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the  Rockwell  C  hardness  measured  on  flats  ground  along  the  side  of  the 
end  quench  test  bar.  This  is  possible  because  there  is  a  pronounced 
hardness  gradient  caused  by  the  microstructural  gradient.  Although  sur¬ 
veying  hardness  is  a  qreat  deal  easier  and  less  time  consuming  than 
directly  observing  the  microstructure  under  a  microscope,  the  hardness 
surveys  are  tedious  because  up  to  60  individual  manual  hardness  measure¬ 
ments  may  be  required  to  survey  a  single  test  bar. 

Figure  1  compares  a  conventional  hardness  survey  with  a  longitudinal 
velocity  scan  of  an  end-quench  test  bar  of  AISI  type  4140,  a  common,  low 
alloy,  heat-treatable  steel.  The  acoustic  path  was  transverse  to  the 
axis  of  the  test  bar.  Parallel  flats  were  ground  on  opposite  sides  the 
full  length  of  the  test  bar;  the  acoustic  path  was  thus  along  bar  diameters 
through  the  thickness  of  the  bar  between  the  ground  flats.  The  relative 
acoustic  velocity  is  plotted  as  Av/vQ  ,  the  fractional  change  in  velocity 
at  any  point  relative  to  vQ  the  velocity  at  the  slower  cooled  pearl itic 
end  of  the  test  bar.  As  seen  in  the  plot,  the  velocity  decreased  by 
about  0.7%  at  the  martensitic,  quenched  end.  The  accuracy  of  the  relative 
acoustic  velocity  measurement  is  limited  by  the  uniformity  of  the  thick¬ 
ness  and  the  deviation  from  perfect  parallelism  of  the  ground  flats. 

This  is  estimated  to  be  1  or  2  parts  in  104  which  is  about  2%  of  the 
actual  range  of  velocities  measured.  Thus  the  sensitivity  and  discrimi¬ 
nation  of  the  velocity  measurement  is  equal  to  or  better  than  that  of  the 
hardness  measurements. 

The  computer  controlled  automatic  acoustic  velocity  scan  requires  only 
4  minutes  from  start  to  finish  with  the  data  automatically  plotted  as 
Av/Vg  versus  position  in  the  bar.  This  contrasts  favorably  with  the 
hour  or  so  required  for  a  manual  hardness  survey.  To  our  knowledge,  this 
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FIG.  l--£nd-quench  test  results  measured  in  type  4140  steel.  The  two  curves  compare  the  micro¬ 
structure  gradient  as  surveyed  acoustically  and  by  conventional  hardness  measurements. 


is  the  first  time  an  end-quench  test  has  been  performed  using  a  velocity 
scan  rather  than  a  hardness  scan. 

In  addition  to  the  4140  steel  end-quench  test,  we  have  also  performed 
similar  acoustic  velocity  scans  on  end-quench  test  bars  of  types  52100, 
4615,  8640,  and  1095  steels  with  similarly  encouraging  results.  We  find 
these  results  persuasive  that  acoustic  velocity  measurements  can  be  used 
to  map  out  microstructural  variations  within  steel  objects. 

C .  _ Acoustic  Attenuation  Measurements 

In  principle,  attenuation  measurements  are  more  attractive  than 
velocity  measurements  as  the  basis  of  acoustic  NDE  microstructural 
characterization.  Where  velocity  is  weakly  a  function  of  microstructure, 
acoustic  attenuation  is  strongly  affected  by  microstructure  because  of 
scattering  at  grain  boundaries,  second  phase  particles,  and  other 
microstructural  features.  Moreover,  the  strong  frequency  dependence  of 
the  acoustic  attenuation  coefficient  can  provide  additional  information 
related  to  microstructure. 

A  sophisticated  technique  has  been  developed  for  automatically 
measuring,  calculating,  and  displaying  the  attenuation  of  acoustic  waves 
through  solid  samples  for  a  wide  range  of  frequencies.  The  technique 
consists  of  launching  a  broadband  pulse  of  acoustic  wave  energy  into 
the  sample  through  a  buffer  and  then  recording  the  echoes  from  the  front 
face,  as  well  as  the  first,  and  the  second  back-face  echo  from  the 
sample.  The  three  pulses  are  then  digitized  and  stored  in  computer 
memory,  and  subsequently  used  as  the  basis  for  introducing  corrections 
due  to  transducer  response,  reflections  from  buffer/sample  interfaces, 
and  for  calculatino  the  attenuation  of  the  material.  The  individual 
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frequency  components  of  the  broadband  pulses  are  corrected  in  amplitude 
for  diffraction  effects,  and  the  attenuation  as  a  function  of  frequency 
is  plotted  direction  by  the  computer.  Signal  averaging  and/or  sampling 
techniques  are  also  used  for  increasing  the  signal-to-noise  ratio,  and 
for  processing  signals  in  the  frequency  range  from  5  to  70  MHz. 

Figure  2  shows  the  attenuation  in  two  different  steel  samples  with 
identical  compositions  but  different  microstructure.  The  curves  for 
pearl ite  and  martensite  show  frequency  dependences  close  to  f^  indi¬ 
cating  that  Rayleigh  scattering  is  the  dominant  mechanism  for  attenuation. 
Pearl itic  data  was  obtained  from  two  different  transducers  as  shown  by 
the  overlap  on  the  figure.  These  first  results  demonstrate  the  ability 
of  our  attenuation  measuring  technique  to  discriminate  between  steel 
microstructures.  A  variety  of  experiments  are  in  progress  to  apply  the 
method  to  steels,  copper,  and  aluminum  alloys.  We  believe  that  attenuation 
measurements  can  be  used  to  gain  much  useful  information  about  the  micro- 
structural  state  of  materials. 


III.  MEASUREMENTS  OF  STRESS  STATES  BY  ACOUSTIC  NDF 


(D.  Barnett,  G.  Herrmann,  W.  r.  J.  Deeg,  and  R.  King) 

A.  Introduction 

In  numerous  aeronautical  and  aerospace  structural  systems  it  is 
imperative  to  have  available  accurate,  reliable  methods  for  periodic 
nondestructive  evaluation  of  critical  parts  in  order  to  predict  useful 
service  lifetimes.  These  nondestructive  inspections  are  needed  to  ensure 
the  integrity  of  structural  elements  under  complex  loading,  thermal,  and 
environmental  conditions,  with  particular  emphasis  on  the  prevention  of 
catastrophic  or  brittle  failure. 

During  the  past  year,  we  have  worked  toward  developing  entirely  new 
techniques  for  quantitatively  determining  the  states  of  stress  associated 
with  cracks  and  other  sources  of  stress  concentration  using  acoustic 
nondestructive  methods.  In  particular,  we  have  made  very  good  progress 
in  developing  methods  for: 

1.  Determining  whether  or  not  stress  concentrations  such  as  flaws. 
Inclusions,  or  cracks  are  present;  and 

2,  Assessing  whether  or  not,  for  given  loading  conditions,  cata¬ 
strophic  failure  due  to  unstable  crack  propagation  is  a  distinct 
possibility,  without  the  need  for  making  direct  measurements  on 
the  stress  concentration  sites. 

At  the  same  time  we  have  recognized  that  under  many  technological  applica¬ 
tions  of  engineering  importance,  real  materials  are  often  subjected  to 
the  effects  of  rather  complex  coupled  fields,  e.g.,  piezoelectric, 
magnetoelastic,  and  thermoelastic  fields.  The  governing  fracture  mechanics 


for  materials  experiencing  coupled  fields  had  not  previously  been  developed, 
and  we  have  attempted,  with  success,  to  construct  the  extended  fracture 
mechanics  framework  required.  We  view  this  as  a  significant  advance, 
since,  just  as  in  the  case  of  purely  elastic  solids,  the  fracture 
mechanics  framework  provides  the  vehicle  for  successful  coupling  of 
theory  and  experiment. 


B. The  Acoustic  Determination  of  the  J  and  M  Integrals  in  Cracked  Plates 

Let  us  consider  an  elastic  solid  in  a  state  of  two-dimensional  plane 
deformation.  The  line  integrals 

9u. 

J  =  $  (Wdy  -  T.  —  ds)  (1) 

c  3x 

9uv 

M  =  j  (Wn.x.  -  TK  — -  x. )  ds  (2) 

c  ax.  ' 

may  be  shown  to  yield  the  generalized  thermodynamic  forces  (so-called 
crack  extension  forces)  on  single-ended  and  double-ended  sites  of  stress 
concentrations  with  the  contour  c  .  Here  W  is  strain  energy  density, 
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T  and  u  are  the  traction  and  displacement  vectors,  respectively,  n 
is  the  outer  normal  to  c  ,  and  ds  and  dy  are,  respectively,  an 
element  of  arc  of  c  at  x  and  its  component  along  the  y-direction. 

If  c  contains  one  tip  of  a  crack  parallel  to  the  x-direction,  then 

2 

1  -  v  p  p 

J  =  — - -  (Kj  +  Kjj)  =  G  (3) 

where  v  and  E  are  Poisson's  ratio  and  Young's  modulus,  respectively. 
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Kj  and  Kjj  are  the  Mode  1  and  Mode  II  crack  stress  intensity  factors, 
and  G  is  the  crack  extension  forces.  Similarly,  one  can  show  that  if 
c  encloses  both  tips  of  a  crack  of  length  £  then 

M  =  £G  .  (4) 

Both  J  and  M  are  path- independent  so  that  the  contour  c  is  quite 
arbitrary  and  may  be  chosen  for  both  theoretical  and  experimental  conven¬ 
ience.  Since  the  usual  fracture  criterion  is  expressed  as 

G  =  Gc  (5) 

where  Gc  is  the  fracture  toughness  which  is  a  material  property  and 
determinable  by  separate  destructive  measurements;  it  is  immediately 
clear  that  if  one  can  determine  J  and  M  nondestructively,  say,  by  acoustic 
methods,  one  can  use  this  information  to; 

1.  Determine  whether  or  not  stress  concentration  sites  exist  in  a 
given  region,  i.e.,  if  both  J  and  M  vanish  there  are  no  such 
sites;  and 

2.  If  stress  concentration  sites  exist,  determine  whether  or  not  G 
is  dangerously  close  to  Gc  . 

During  the  past  year  we  have  shown  that  it  is  indeed  possible  and 
practical  to  determine  J  and  M  from  acoustic  measurements  on  loaded 
cracked  panels  from  measurements  of  the  relative  velocity  shift  of 
longitudinal  and  shear  waves  propagating  through  the  panels.  It  is  known 
that  longitudinal  and  shear  wave  scans  of  stressed  samples  allow  one  to 
determine  the  in-plane  principal  stresses  Oj  and  o£  as  well  as  the 
orientation  0  of  the  principal  stress  axes  at  selected  points  in  such 
samples.  During  the  past  year  we  performed  an  analysis  (R.  King,  G.  Kino, 
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and  D.  Barnett,  in  preparation)  which  shows  that  J  may  be  expressed  as 
1  -  v2  , 

J  = - p  (o^  .  a  )(a  +  o^)  cos  (20  -  <f>)  ds 

2E  c 

-  1/2  $  o)((oj  -  o^)  sin  (24>  -  G)  +  (oj  +  o^)  sin  0)  ds  .  (6) 

All  quantities  in  the  integrands  in  (6)  are  determinable  directly  from 
acoustic  data  except  for  the  rotation  field  w  ,  which  can  be  found  from 
forward  integration  of  the  compatibility  equations  expressed  in  terms 
of  the  acoustically  measured  stress  state.  A  result  similar  to  (6) 
holds  for  the  M  integral.  Thus,  in  principle,  both  J  and  M  may  be  found 
from  nondestructive  acoustic  measurements  alone. 

For  the  past  two  years  we  have  had  and  are  refining  a  scanning 
longitudinal  wave  system  capable  of  performing  the  requisite  acoustic 
measurements;  development  of  a  scanning  shear  wave  system  is  now  well 
underway.  During  the  past  year  we  have  demonstrated  by  computer  simula¬ 
tion  using  exact  elastic  fields  with  superimposed  random  noise  to  simulate 
experimental  error  to  show  that  acoustic  determination  of  J  and  M  is 
indeed  feasible. 

It  has  also  been  clear  that  a  successful  program  to  acoustically 
determine  J  and  M  must  involve  comparison  with  theory.  To  this  end, 
one  requires  theoretical  knowledge  of  the  exact  elastic  fields  at  points 
in  geometrical ly  complex  plane  laboratory  samples  under  external  loading. 
Such  complete  solutions  are  almost  never  available  in  the  literature. 
During  the  past  year  we  have  developed  a  Boundary  Integral  Equation 
computer  code  capable  of  solving  problems  of  cracked  elastic  plates  of 
arbitrary  geometry  and  of  arbitrary  elastic  anistropy.  The  technique 
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will  not  be  described  here,  but  we  note  that  development  of  such  codes 
has  been  supported  in  other  laboratories  by  other  branches  of  AFOSR. 

To  our  knowledge,  ours  is  the  first  Boundary  Integral  Code  valjd  for 
arbitrary  anistropy.  The  computer  program  is  quite  compact  and  user- 
oriented.  A  paper  describing  the  program  (D.  Barnett  and  R.  King)  is 
currently  being  prepared. 

C.  Stress  Analysis  and  Fracture  Mechanics  in  Solids  Containing  Cracks  and 
Inclusions  and  Subjected  to  Coupled  Fields 

We  have  generated  the  first  complete  solution  for  the  stress  analysis 
and  fracture  mechanics  of  a  slit-like  crack  in  a  piezoelectric  solid  of 
arbitrary  anistropy  (piezoelectrics  are  necessary  anistropy).  This  work 
is  currently  being  prepared  for  publication  (W.  F.  J.  Deeg  and  D.  Barnett, 
to  be  submitted  to  International  Journal  of  Sol  ids  and  Structures) .  The 
stress  analysis  developed  shows  how  both  externally  applied  loads  and 
electric  fields  influence  the  stress  intensity  factors  and  the  crack 
extension  force  in  such  media.  A  remarkable  result  of  the  analysis  is 
that  the  crack  extension  force  may  be  either  positive,  negative,  or  zero 
depending  upon  the  ratio  of  applied  tractions  to  the  applied  electric 
field!  In  essence,  this  means  that  given  a  cracked  piezoelectric  with 
certain  applied  loadings,  the  crack  may  be  rendered  either  stable  or 
unstable  by  varying  the  applied  electric  field  appropriately.  The  exten¬ 
sion  of  the  J  and  M  integrals  to  such  solids  has  also  been  made. 

An  analysis  with  even  more  far-reaching  consequences  has  been  developed 
for  problems  of  piezoelectrics  containing  voids  and  inhomogeneities  has 
been  completed  (W.  F.  J.  Deeg  and  D.  Barnett,  to  be  submitted  to  Inter¬ 
national  Journal  of  Engineering  Science).  We  have  extended  a  previous 
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analysis  by  J.  0.  Eshelby  for  "transformed"  elastic  inclusions  and 
inhomogeneities  to  the  piezoelectric  case.  These  are  fully  three-dimen¬ 
sional  anistropic  solutions.  The  problems  concerned  the  fields  associated 
with  inclusions  and  inhomogeneities  which,  in  the  absence  of  a  surrounding 
matrix,  would  like  to  take  on  a  stress-free  or  "transformation"  strain. 

The  constraintof  the  surroundings  induces  an  initial  stress  and  strain 
state  in  the  entire  medium.  The  transformation  strain  may  arise  from  a 
phase  transformation  or  be  due  to  thermal  mismatch;  a  slight  re-interpreta¬ 
tion  of  the  transformation  strain  allows  us  to  completely  treat  the  problem 
of  stress  concentrations  about  ellipsoidal  voids  and  inclusions  perturbing 
far-field  applied  loadings  and  electric  fields  in  a  very  simple  way. 

Hence,  we  now  feel  that  we  have  made  a  rather  complete  development  for 
treating  problems  of  stress  concentrations  in  infinite  piezoelectric  solids. 

For  piezoelectric  solids  of  finite  dimensions  and  complex  boundary 
shapes  we  have  formulated  the  appropriate  Boundary  Integral  Equation 
method  for  generating  the  solutions  of  interest  to  us.  This  involves 
a  minor  overhaul  of  our  previously-mentioned  Boundary  Integral  Code  for 
pure  elastic  solids,  and  we  anticipate  that  the  computational  scheme  will 
be  fully  operational  within  the  coming  year. 

Finally,  we  note  that  the  notion  of  conservation  integrals,  the  0  and 
M  integrals,  for  thermoelastic  and  porous  solids  has  been  treated  during 
the  last  year. 
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Abstract 

Based  on  invariant  variational  principles  due  to  Noether,  it  is  possible 
to  derive  several  conservation  theorems  of  mathematical  physics  applicable, 
in  principle,  for  any  field  theory.  These  conservation  laws,  and  the  resulting 
so-called  path- independent  integrals,  have  proved  useful  in  continuum  mechanics 
and  in  particular  in  the  theory  of  elasticity. 

Several  new  applications  of  these  conservation  theorems  will  be  discussed, 
including  specifically  those  to  elastodynamics ,  thermoelasticity,  porous  media, 
and  the  experimental  determination  of  stress  intensity  factors  in  fracture 


mechanics. 


1. 


INTRODUCTION 


A  celebrated  theorem  due  to  Emmy  Noether  (1882-1933)  [ 1 J  makes  it  pos¬ 
sible  to  construct  divergence-free  quantities  from  the  variables  of  a  varia¬ 
tional  problem,  provided  the  extremum  integral  possesses  certain  invariance 
properties. 

If  Hamilton's  Principle  is  applied  to  derive  varint ionally  the  equations 
of  motion  of  a  physical  theory,  one  obtains  merely  the  Euler-Lagrange  equations 
because  the  domain  ("fixed  endpoint  variation")  itself  is  not  varied. 

By  contrast,  if  the  same  variational  integral  (whose  integrand  is  the 
Ldgrangian  density)  is  subjected  to  a  more  general  type  of  variation,  allow¬ 
ing  endpoint  variations,  certain  additional  terms  do  appear,  along  with  the 
terms  corresponding  to  the  Euler-Lagrange  equations. 

The  additional  terms  form  a  tensor,  referred  to  as  the  energy-momentum 
tensor.  Its  usefulness  in  continuum  mechanics,  and  particularly  in  the  theory 
of  elasticity,  as  discussed  by  Eshelby  in  a  series  of  papers  beginning  in 
1951  [2,  3,  4],  rests  on  the  fact  that  the  energy-momentum  tensor  leads  to 
conservation  integrals  (so-called  "path-independent"  integrals)  which  have 
an  irmediate  physical  significance.  Limited  availability  of  space  does  not 
pe~.i:  here  a  detailed  discussion  of  these  integrals.  It  will  suffice  to  men¬ 
tion  that  the  so-called  translation  integral,  first  discussed  by  Eshelby  [2], 
represents  a  force  on  a  defect.  A  defect  is  defined  in  this  context  as  any 
type  of  nonhomogeneity,  inclusion,  cavity,  dislocation,  crack,  etc.  The 
"force"  on  a  defect  is  not  the  customary  Newtonian  force,  but  is  defined, 
rather,  as  the  rate  of  change  of  total  energy  of  the  system  with  respect  to  a 
unit  infinitesimal  displacement  of  the  defect. 

Quite  independently  of  Eshelby,  three  conservation  laws  of  clastostatics 
were  derived  by  Gunther  15).  lie  showed  how  all  three  of  them  could  be  applied 
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in  problems  of  torsion,  plane  stress,  beam  theory,  plate  theory,  and  shell 
theory.  Gunther  was  not  interested  in  defects,  but  simply  in  providing  ex¬ 
pressions  of  boundary  integrals  which  were  path- independent .  It  can  be  rig¬ 
orously  shown  that  in  Euclidean  space  there  will  be  just  three  conservation 
laws,  because  only  three  groups  of  infinitesimal  transformations  satisfy 
Noether's  invariance  principle.  The  three  conservation  laws  are  associated 
with  che  transformation  of  translation,  rotation  and  similarity  (or  gauge), 
respectively.  The  translation  integral,  which  provides  Eshelby's  force  on 
a  defect,  was  independently  deduced  by  Rice  [6]  for  the  special  case  of  a 
crack,  which  he  called  the  J-integral.  The  rotation  and  similarity  integrals 
were  independently  discovered  by  Knowles  and  Sternberg  [7],  who  did  not  dis¬ 
cuss  their  potential  application,  but  extended  them  to  nonlinear  elasticity. 

In  the  context  of  fracture  mechanics,  the  importance  of  path- independent 
integrals  resides  in  the  fact  that  they  can  be  related  to  energy  release 
rates,  i.e.  to  crack  extension  forces  which  themselves  depend  only  on  stress 
intensity  factors. 

It  was  shown  by  Rice  [6]  that  the  J-integral  is  related  to  the  total 
potential  energy-release  rate  associated  with  moving  or  extending  cracks  (in 
their  ova  plane)  in  linear  or  nonlinear  elastic  materials.  Later,  Budiansky 
and  Kice  (S)  have  shown  that  the  two  additional  integrals  of  Gunther  and  of 
Knowles  and  Sternberg  are  associated  with  cavity  or  crack  rotation  (they  called 
it  the  L-integral)  and  expansion  (or  similarity,  they  called  it  the  M-integral), 
respectively . 

In  what  follows,  the  energy-momentum  tensor  is  briefly  derived  and  applied 
to:  a)  free  motions  of  an  clastic  bar;  b)  equilibrium  of  a  thermoplastic  solid 
or  fluid-filled  porous  solid  and  c)  to  the  experimental  determination  of  stress- 
intensity  factors,  making  use  of  the  J-integral. 
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ENERGY-MOMENTUM  TENSOR 


As  suggested  by  Eshelby  [3, A)  and  by  Haul  [9),  the  energy-momentum  tensor 
can  be  derived  in  a  concise  fashion  without  explicit  recourse  to  the  calculus 
of  variations.  Let  the  Lagrnngian  density  L  depend  on  n  independent  var¬ 
iables  x°  ,  on  m  dependent  functions  <>*(xa)  ,  as  well  as  on  the  gradients 
id  i 

3$  /?:.  E  .  The  Lagrangian  density  can  also  be  considered  a  function 

of  only  the  independent  variables,  which  will  be  denoted  by  L  .  Thus  we 
can  write 

L  =  L(x°  ,  ,  3  $*)  *  £(xQ)  a  *=  1,2,. . .  ,n  ;  i  =  1,2,. . .  ,m  (2.1) 

a 

If  we  define 

3L/D*1  E  Ti  ;  SL/a^1  E  T.a  (2.2);  (2.3) 

then  the  chain  rule  will  give 

-  >L  +  v/  +  <2-« 

which  can  also  be  written  as 


3L'3L  +  T.?61+  3  (T  03  ^.1)  -  (3  T  P)  3  ^ 
a  a  i  a  Baa  3  i  a 


(2.5) 


or 


■  9aL  +  (Ti  -  Vi”’  V*  +  VTiV> 


(2.6) 


The  parenthesis  T^  -  3^T^  is  recognized  to  be  the  i-th  Lagrange  equation 
of  motion  which  will  be  abbreviated,  ns  usual  by  IL]^  .  If  the  last  term  in 
eqn.(2.6)  is  brought  to  the  left-hand  side,  this  equation  may  be  rewritten  at 


vV  -  ■  v-*  ^i9/ 


(2.7) 
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The  negative  of  the  quantity  in  parenthesis  is  defined  as  the  energy-momentum 

O 

tensor  P 

a 

P  B  =  -£  6  &  +  T,B3  (f.1  (2.8) 

a  a  i  a 

If  all  the  independent  variables  are  cyclic  (or  iguorable) ,  i.e.  J  L  =  0  , 
and  if  the  equations  of  motion  are  satisfied,  i.e.  [ LJ  ^  =  0  ,  then  we  see  that 

*0P  3  -  0  (2.9) 

8  a 

i.e.  the  energy-momentum  tensor  is  divergence-free.  This  implies  the  exis¬ 
tence  of  a  potential 


„  YP  &P-,  .  c 

e  e  3  <s 
a  o  YP  P 


(2.10) 


It  may  be  remarked  further  that  the  energy-momentum  tensor  is  the  Hamiltonian 
of  classical  dynamics,  where  one  is  concerned  with  only  one  independent  var¬ 
iable  (time).  Similarly  as  in  classical  dynamics,  one  can  derive  in  the 
present  case  of  n  independent  variables  the  canonic  equations  which  has 
been  done  by  Kane  [10]. 

3.  ELASTIC  BAR 

As  a  simple,  but  not  trivial  example,  let  us  consider  a  continuum  which 
would  involve  two  independent  variables  (and  thus  be  essentially  different 
from  classical  dynamics)  and  only  one  dependent  variable.  One-dimensional 
motions  in  an  elastic  homogeneous  bar  represent  such  a  system.  The  Lagrangian 
density  in  usual  notation  is  given  by 


»  O  9 

2L  =  pu*  -  Eu’ 


(3.1) 


Here  p  is  the  mass  density  per  unit  of  length  of  the  bar,  E  is  Young's 
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modulus  of  the  material,  ti  is  the  displacement.  Dot  am!  prime  Indicate 
partial  differentiation  with  respect  to  time  and  the  x-coordlnato  along  the 
bar,  respectively. 

In  terms  of  the  notations  introduced  in  the  previous  sect  ion,  wo  cm)  write 
i  *  1  ;  a  »  1,2  ;  -  x  ;  x^  »  t  ;  u  ■  ;  u*  **  3^'  ;  u  -  32^* 

The  components  of  the  energy-momentum  tensor  in  this  example  are  calculated 
to  be 

P^  >'  ou*/2  +  pu/2  ,  l'j^  *  -pu*  ,  Pj,^  "  ou  ,  P  ^  “  -ou*/2  —  pu/2  (3.2) 

where  o  is  the  stress  and  p  the  momentum 

o  -  Eu’  ;  p  “  pu  (3.3) 

Since  L  does  not  depend  on  x  and  t:  explicitely,  and  if  the  equation  of 
motion  is  satisfied,  the  energy-momentum  tensor  is  divergence-free  which  leads 
in  this  example  to  the  two  Noether  equations,  in  usual  notation 

3  (ou’/2  +  pu/2)  -  3  (pu*)  »  0  (3.4) 

A  l 

9(ou)  ~  3  (ou’/2  +  pu/2)  “  0  (3.5) 

It  is  tempting  to  try  to  provide  a  physical  int reprotat ion  to  the  above  two 
equations.  The  second  equation  is  more  familiar  in  this  regard.  The  produet 
ou  is  recognized  to  represent,  the  energy  flux  across  a  cross-section  of  the 
bar.  Since  the  sum  ou*/2  +  pu/2  represents  the  total  energy  density  of  the 
bar,  the  second  equation  states  that  the  time-rate  of  change  of  the  total 
energy  of  a  bar  element  has  to  equal  the  not  energy  flux  within  the  element  . 

Thus  equation  (3.5)  Is  a  statement  of  the  principle  of  energy  conservation. 

To  assign  a  physical  meaning  to  equation  (3.4)  Is  not  so  straight-forward. 
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If  it  is  recalled  that  the  force  on  a  defect  is  defined  as  the  spatial  rate 
of  change  of  the  total  energy,  the  first  term  can  he  looked  upon  as  a  kind 
of  force,  which  could  be  termed  non-Newtonian  or  quasi-force.  The  second 
term  has  been  referred  to  in  the  past  ( 3 J  as  quasi-momentum,  pseudomomontum, 
or  field  momentum.  Thus  equation  (3./*)  appears  to  be  a  statement  of  the 
principle  of  linear  momentum,  however  not  in  the  ordinary,  Newtonian  sense, 
but  rather  pertaining  to  the  motion  of  defects  with  respect  to  the  body, 
for  a  further  discussion  see  Eshelby  (3,4]  and  Kogula  (11). 

In  symbolic  notation  equations  (3.4)  and  (3.5)  can  be  written  compactly  as 

V  *  IP  -  0  where  ^  5  ^  *[x  +  (3.6) ;  (3.7) 

and  |P  is  the  energy-momentum  dyadic  with  components  given  in  equation  (12). 
The  divergence  theorem  loads  from  equation  (3.6)  immediately  to  equation 

$  |n  •  \T  ds  =>0  (3.8) 

where  ds  is  the  element  of  length  of  any  closed  circuit  in  the  x,t  plane 
and  |n  is  the  outward  drawn  normal  unit  vector.  As  a  simple  application 
of  the  path- independent  integral  (3.8)  let  us  consider  a  bar  of  length  £ 
and  a  rectangular  circuit  as  indicated  in  the  sketch,  where  t*  is  some 
arbitrary  time.  Equation  (3.8)  leads  then  to  the  relation 

i  1  t*  t* 

/oudx  |  «  Jo udx  |  -  Jpu'dt  |  -  JpuMt  |  (3.9) 

o  tat*  o  t^o  o  x“£.  o  X“o 

The  first  term  on  the  right-hand  side  of  equation  (3.9)  can  ho  evaluated 
on  the  basis  of  prescribed  initial  conditions.  Similarily,  the  second  and 
third  term  on  the  right-hand  side  can  be  evaluated  on  the  basis  of  prescribed 
boundary  conditions.  Relation  (3.9)  may  be  readily  interpreted  as  stating 
that  the  total  (in  spare)  energy  flux  (or  power)  is  conserved  in  time. 
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4.  THERMO ELASTIC  SOLID  OR  FLUID-FILLED  POROUS  SOLID 

Variational  principles  in  thermoo.last  icily  and  heat  conduction  were  dis¬ 
cussed  by  Herrmann  in  112}.  Using  the  analogy  established  by  lUot  (13],  the 
following  development  is  valid  also  for  an  elastic  porous  solid,  whose  pores 
are  filled  with  a  compressible  viscous  fluid. 

The  equilibrium  equations  can  be  written  down,  in  the  absence  of  body 
forces  as 


’y,j  -  V.j  ' 0 


°.i  -  “i  "  0 


and  the  constitutive  equations  as 


Y..-W  =0  ;  Y  +  Wn-’C)  ;  s  +  D  -  0  (4.3) ; (4.4) ; (4,5) 

In  these  equations  t  -  P^G  ~  °i }  denotes  the  total  stress,  t^j  is  the 

isothermal  part  of  stress,  P^j  is  related  to  the  thermal  expansion  properties 

of  the  material,  G  is  the  temperature  increment  above  a  reference  absolute 

temperature  T  ,  g  is  the  thermal  disequilibrium  force,  y  e  (u  + 

1  l  ij  i *j 

u,  )/2  is  the  clastic  strain,  yb  s.  ,  +  P,,u  ,  is  the  "thermoelastic 
dilatation".  The  entropy  displacement  s_j  is  related  to  the  entropy  s 
by  the  equation 


s  -  -s1#i  (4.6) 

W  is  Biot's  thcrmoelastic  potential,  expressed  here  as  a  quadratic  form  of 
isothermal  components  of  stress  and  the  temperature  increment 

2W  "  BiJkfcTijTk*.  +  C° 

whore  c  Jo  tlio  specific  heat  per  unit  volume  for  zero  strain.  D  is  Biot's 
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dissipation  function,  expressed  here  us  a  quadratic  form  of  the  "thermal  dis¬ 
equilibrium  force"  g^  ,  the  quantity  conjugate  to  the  entropy  displacement  s^. 


“  ‘  klJ*i*j/pTr 


0.8) 


Here  is  the  thermal  conductivity  tensor  and  p  is  the  time  operator 

3/3t  .  To  derive  a  variational  principle  for  displacements,  the  const itutuve 
equations  are  written  in  inverted  form 


I,,  c  U 

ij  ,Y 


-U 


1J 


>Y 


Ei  -  G 
i  ,Si 


(A.9) ;  (4.10) ; (4.11) 


where  the  quadratic  forms  U  and  C  are 


U  “  2  cijkiYijYk.».  +  2  TrY  /c 

r  iPTy 
g  “  2  ir~  sisi 
ij  1 


(4 • 12) 


(4.13) 


Consider  now  the  Lagrangian  density  L  to  be  U  -  G  and  let  the  independent 


variables  be  y..  ,  y  and  s.  ,  i.e. 

ij  i  ’ 


L  =  U  -  C.  *  My^  ,  y  ,  sj.)  (4.14) 

The  thermoelastic  energy  momentum  tensor  can  then  be  calculated  and  can  be 
represented  in  either  one  of  the  following  forms 


Pij 

(U  -  G  +  Oy)  6^  -  Tkjuk>1 

(4.15) 

Pij 

^2  akiYkt  +  2  Y°+  2  Rkr,k^  ^ij  “  TkjUk,i 

(4.16) 

PU 

(U  "  C  "  °S)  *ij  "  °kj\,i 

(4.17) 

where  o.  ,  «  x.  ,  -  p.  ,0 
kj  kj  kj 

(4.18) 

Since  the  energy-momentum  tensor  is  divergence-free,  a  "path-independent" 
integral  can  immediately  be  constructed  formally.  In  the  absence  of  thermal 
effects  it  will  reduce  to  Rice's  J-integral.  We  obtain  either  one  of  the 
following  forms 


Jk  -  /  l  (U  -  G  - 
Jk  "  /  l2(xijYij 
Jk  "  /  l2(aijYij 


s)  nk  -  n^uj 

Y°  +  Vi*  nk 

S0  +  s  g  )  nk 


k1  da 

Vijuj,k]  da 

T j  ° j  ,  k^  dn 


(4.19) 

(4.20) 

(4.21) 


Integration  is  to  be  carried  out  over  a  closed  surface  with  element  da  . 
Wherever  the.  time  operator  p  =  -r—  occurs,  the  multiplication  has  to  be 

O  u 

carried  out  as  convolution,  similarly  as  in  the  dynamic  case  discussed  by 
Gurtin  (14]. 

5.  EXPERIMENTAL  DETERMINATION  OF  STRESS-INTENSITY  FACTORS 

In  some  recent  and  still  ongoing  experimental  studies,  the  stress-inten¬ 
sity  factor  at  the  root  of  an  edge-notched  6061-T6  aluminum  specimen  in  ten¬ 
sion  has  been  determined  by  direct  measurement  of  the  J-integral,  defined  in 
two  dimensions  as 

3u 

J  "  J  (Wdy  -  Tt  ds)  (5.1) 

Here  T  is  a  curve  surrounding  the  crack  tip,  W  the  energy  density,  x  and 
y  are  rectangular  coordinates  parallel  and  normal  to  the  crack,  and  ds  is 
an  clement  of  arc  length  along  I’  ,  T^  is  the.  traction  vector  defined  with 
respect  to  the  outward  normal  along  I’  ,  uA  is  the  displacement  vector. 

The  value  of  J  has  been  determined  experimentally  by  measuring  W  ,  T^ 
and  Du^/dx  on  the  rectangular  contour  of  the  test,  specimen.  The  specimen 
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was  1.0  cm  thick,  6.0  cm  wide  and  the  distance  between  shoulders  was  15.2  cm. 
The  edge  crack  was  placed  midway  between  shoulders  and  had  a  depth  of  0.96  cm. 
If  the  specimen  is  thought  to  be  in  vertical  position,  then  along  the  ver¬ 
tical  edges  0  and  the  contribution  to  J  is  only  from  f  Wdy  .  Across 

the  top  and  bottom  of  the  sample  dy  =  0  and  ds  *=■  dx  such  that  only 
-/  T^(?u^/?x)  dx  contributes  to  J  .  To  evaluate  W  ,  10  strain  gages  were 
placed  along  the  upper  half  only  (due  to  symmetry)  of  the  vertical  edges  and 

it  was  assumed  that  the  specimen  was  linearly  elastic  everywhere,  thus 
2 

W  »  Ec  /2  where  E  is  Young’s  modulus  and  e  is  the  strain.  Further,  W 
is  proportional  to  P  ,  the  load  on  the  sample.  The  dimensions  of  the  speci¬ 
men  were  chosen  such  that  the  strain  just,  below  the  shoulders  of  the  tension 
specimen  would  he  approximately  uniform.  This  au^/^x  could  be  determined 
with  good  accuracy  using  only  several  points  across  the  specimen,  obtaining 
nearly  the  same  value  (rigid  rotation). 

The  contribution  to  J  from  the  vertical  edges  (strain  gage  measure¬ 
ments)  was  determined  to  be 

J  -  0.252  x  10"l7P2E  (5.2) 

N 

The  contribution  to  J  from  the  horizontal  portion  of  the  contour,  measured 
by  means  of  linear  variable  differential  transformers  (LVDT)  as  displacement 
gages  was  found  to  be 

J  »  1.96  x  10~17P2E  (5.3) 

IT 

Addition  gives  the  total  contribution  to  .1 


J  ”  JV  +  JH  "  2,21  X  -10~17,’2e  ^ 


(5./.) 
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2  j 

The  equivalent  stress  intensity  value  is  K  •*  [JE/(1  -  v  ))’  for  plane 
strain  cracks.  Thus  the  experimental  method  permitted  to  measure  K/F,  . 

This  can  be  compared  with  the  calculated  handbook  value  115),  which  is 

i  ~  , ,  , „-17„2  m 

J  »  2.15  *  10  P  e  — =• 

N 

10  2 

Here  a  value  of  E  ■  7.03  *  10  N/m  has  been  used. 

It  is  seen  that  the  agreement  between  the  theoretical  and  the  measured 
value  of  J  is  excellent..  Note  that  only  about  10%  of  the  value  of  .1 
stems  from  the  vertical  part  of  the  integral  for  the  specimen  with  the  de¬ 
scribed  geometry.  A  more  detailed  account  of  this  invest  igat  ton  is  given 
in  ref.  [16]. 
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Abstract 

We  arc  using  measurements  of  acoustic  velocity 
and  attenuation  for  nondestructi ve  characterization 
of  the  microstructure  of  materials.  Ihe  spatial 
variation  of  the  acoustic  velocity  in  Jominy  end- 
quench  spec  fens  was  found  to  correlate  very  well 
with  Rockwell  C  hardness  scan:,  ir.aicating  a  poten¬ 
tially  practical  method  for  measuring  the  hardening 
response  of  heat  treated  steel.  Attenuation  meas¬ 
urements  have  also  been  performed  by  using  broad¬ 
band  acoustic  pulses  corrected  for  transducer  res¬ 
ponse,  liquid  buffer/solid  specimen  reflections, 
and  diffraction  effects.  Higher  attenuation  was 
observed  for  pearl itic  than  for  martensitic  micro- 
structures. 


Mechanical  properties  of  materials,  such  as 
tensile  strength,  hardness,  and  impact  toughness, 
are  strongly  dependent  on  microstructural  features 
liLe  grain  size  and  shape,  the  proportions  and 
spatial  distribution  of  the  phases  present,  and 
macroscale  heterogeneities  inmicrostructure  result¬ 
ing  from  local  compositional  or  thermal-mechanical 
processing  variations.  This  is  true  of  both  metal¬ 
lic  and  ceramic  materials.  The  microstructural 
character  of  materials  and  their  related  properties 
are  usually  assessed  by  microscopic  examination  or 
by  direct  measurement  of  properties;  such  proce¬ 
dures  arc  often  costly,  and  are  necessarily  des¬ 
tructive,  requiring  sacrificial  example  specimens. 

Steels  are  alloys  of  iron  and  carbon  modified 
by  minor  additions  of  other  elements,  which  exhibit 
remarkable  variability  of  mechanical  properties. 
These  useful  properties  are  dependent  upon  the 
rather  complex  microstructure  of  steel  which  can 
be  closely  controlled  by  proper  selection  of  com¬ 
position  and  heat  treatment. 

The  dependence  of  steel  microstructure  on  com¬ 
position  and  heat  treatment  is  illustrated  in  Fiq. 
1.  figures  1(a)  and  1(b)  compare  the  microscopic 
appearance  of  polished  and  etched  stool  specimens 
containing  0.?  and  0.6  wt.  pet.  carbon  respect¬ 
ively;  after  annealing  at  900"  C  and  slow  cooling 
in  air  to  room  temperature  both  steels  consist  of 
ferrite,  the  white  grains,  and  pearl ite,  the 
darker  constituent.  Icrritc  is  essentially  pure 


iron;  poarlitc  is  a  mixture  of  two  different  crys¬ 
talline  phases  (too  finely  dispersed  to  be  resolved 
at  the  magnification  shown),  iron  carbide  (IcgC) 
and  ferrite.  Since  virtually  all  the  carbon  res¬ 
ides  in  the  pearl ite,  the  proportion  of  pearl ite 
is  greater  in  the  higher  carbon  steel,  figure  1(c) 
shows  a  steel  with  0.1  wt.  pet.  carbon  with  a 
martensic  microstructuro  rather  than  pearl ite  and 
ferrite.  This  steel  lias  been  cooled  rapidly  in 
water  after  annealing  at  900*0;  the  fast  cooling 
causes  the  formation  of  martensite,  a  supersatur¬ 
ated  solid  solution  of  carbon  in  iron.  Mechanical 
properties  arc  strongly  affected  by  the  proportions 
of  ferrite  and  pearl itc  and  by  the  existence  of 
martensite  rather  than  pearl itc  and  ferrite. 

Changes-  in  microstructuro  such  as  these,  enable 
steels  to  be  heat  treated  to  obtain  optimum  com¬ 
binations  of  properties. 

Our  measurements  of  longitudinal  acoustic  wave 
velocity  in  metallic  samples  \  -  re  evade  with  a  com¬ 
puter-controlled  system  deve!  -d  for  measuring 
acoustic  velocity  fields  in  .!  samples  immersed 
in  a  liquid  buffer  and  descrv  in  detail  in  an 
accompanying  paper  at  this  cc.  rence.* 

1.  Absolute  Velocity  Measurements 

Specimens  of  5  plain  carbon  steels  were  pre¬ 
pared  by  machining  1  cm  thick  flat  plates  oriented 
both  longitudinally  and  transversely  relative  to 
the  steel  bar  stock.  Since  the  acoustic  wave  prop¬ 
agation  direction  was  through  the  1  cm  plate  thick¬ 
ness,  acoustic  velocity  could  bo  measured  both 
parallel  to  and  transverse  tr  the  rolling  direc¬ 
tion  for  each  steel  composition  to  detect  the 
influence  of  any  preferred  orientation  or  crystal 
texturing.  The  compositions  of  the  steels  are 
shown  in  Table  1.  All  specimens  were  given  iden¬ 
tical  heat  treatment;  they  were  heated  to  1<00"C 
and  air  cooled  (normalized),  lhis  resulted  in 
pearl itc/ferritc  microstructures,  the  pearl ite 
varying  front  1?  to  100  vlumc  pet.  over  the  range 
of  carbon  contents.  After  machining,  the  flat 
specimens  were  lapped  to  assure  that  their  flat 
surfaces  were  parallel  to  within  ?.!>  pm. 

The  resulting  velocity  measurements,  which 
include  correction  for  diffraction  effects, 
exhibited  random  scatter  that  totally  obscurred 
the  variation  of  velocity  anticipated  Item  the 
variation  in  microstructuro.  Pig.  2.  One  can 
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easily  calculate  an  expected  acoustic  velocity  from  the 
known  densities  of  pure  Iron  ferrite*  and  fejC,3  and 
known  elastic  constant  data  for  ferrlte/lcjC  mixtures.’ 
Increasing  carbon  should  cause  a  decrease  In  the  longi¬ 
tudinal  acoustic  velocity  linearly  proportional  to 
carton  content;  velocity  should  decrease  0.U6  pet.  per 
wt.  pet.  carbon  in  the  steel,  all  other  factors  being 
Identical.  The  scatter  observed  in  the  mcasurrd  vel¬ 
ocity  was  not  caused  by  variations  In  preferred  orien¬ 
tation;  had  there  been  any  significant  degree  of  pre¬ 
ferred  orientation,  the  longitudinal  and  transverse 
velocities  would  differ  much  wore  than  the  slight 
variations,  seen  in  fig.  2.  We  believe  that  the  ran¬ 
dom  variations  in  such  elements  as  manganese,  silicon, 
sulfur,  and  others  always  present  in  at  least  trace 
amounts,  cause  variations  in  density  and  elastic  mod¬ 
ulus  and  thus  introduce  random  variations  in  acoustic 
velocity  equal  to  or  greater  than  the  systematic 
effect  of  carbon  in  changing  the  nicrostructurc.  For 
cxa.-olc,  \;c  have  calculated  that  the  range  of  manganese 
contents  in  r ur  specimens  (see  fable  I)  causes  density 
variations  large  enough  to  obscure  all  the  above  change 
attributable  to  variations  in  microstructurc.  There¬ 
fore  we  conclude  chat  absolute  velocity  measurements 
arc  not  practical  means  for  characterizing  steel  micro¬ 
structure,  because  random  compositional  variations  will 
confuse  any  attempt  to  calibrate  absolute  acoustic 
velocity  variation  with  wicrostructure. 


(C) 


Fig.  1.  Different  steel  microstructurcs  resulting  from 
difference  in  composition  and  heat  treatment, 
(a)  pearl i to  and  ferrite  in  normalized  0.2  wt. 
pet.  C  steel;  (b)  pearl ite  and  ferrite  in  nor¬ 
malized  0.6  wt.  pet.  C  steel;  (c )  martensite 
in  water  quenched  0.4  wt.  pet.  0  steel 

TABLE  1 


Composition  of  Plain  Carbon  Steels  used  for 
Acoustic  Velocity  Measurements 


Steel  Type 
(AT  SI ) 

C 

Hn 

Si 

P 

s 

1010 

0.10 

0.45 

0.03 

0. 007 

0.032 

•  1020 

0.20 

0.52 

0.21 

0.009 

0.028 

1035 

0.33 

0.72 

0.17 

0.018 

0.025 

1060 

0.55 

0.80 

0.19 

0.012 

0.023 

1095 

0.93 

0.50 

0.22 

0.009 

0.030 

MEASURED  ACOUSTIC  VELOCITY 

•  *  TRANSVERSE  TO  ROLLING 

*  s  PARALLEL  TO  ROLLING 


c> 

V) 

E 

u 


A 


■ 


A 


m 

o 


>- 

H 

o 

o 

to 

> 


5.04 


5.9 
0.00 


0.2.0  0.40  0.60  0.60 

CARBON  (PERCENT ) 


1.00 


Fig.  2.  Measured  acoustic  velocity  in  plain  carbon 

Steels  of  varying  carbon  content.  Tbo  random 
scatter  caused  by  minor  fluctuations  in  alloy¬ 
ing  elements  obscures  the  systematic  variation 
In  velocity  related  to  the  microstructure. 
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t.  Relative  Velocity  Measurements 

fiy  relative  velocity  measurements  we  mean 
measurements  of  acoustic  velocity  at  different 
locations  within  a  given  object.  If  the  material 
cev»t>osition  is  uniform  within  the  object  (as  is 
usually  the  case),  variations  in  acoustic  velocity 
can  be  unequivocal ly  related  to  microslructural 
variations  frooi  one  location  to  another.  There 
are  many  situations  where  microstructure  is  not 
uniform  and  the  pattern  of  microstructural  vari¬ 
ation  is  irportant,  c.g.,  in  an  induction  hardened 
Steel  shaft  teat  treated  to  be  martensitic  on  the 
outer  surface  but  pearl itic  in  its  interior. 

Jominy  end-quench  test  bars  provide  another  example 
of  steel  objects  containing  microstructural  grad¬ 
ients.  We  successfully  used  relative  velocity 
measurements  to  survey  quantitatively  the  micro¬ 
structure  of  cnd-quench  test  bars  of  several  dif¬ 
ferent  alloy  steel  compositions. 

The  Joniny  Test  or  Fnd-Quonch  Test  is  the 
stanJard  metallurgical  quality  control  procedure 
for  measuring  the  hardening  response  of  heat  treat¬ 
able  steels.  Ihc  cnd-quench  test  specimen  is  a 
one- inch  diameter  round  bar.  four  inches  long.  The 
bar  is  heated  to  about  8SC°C.  It  is  then  placed  in 
a  fixture  ae.d  cooled  by  a  jet  of  cold  water 
impinged  upon  one  end  as  sl.oi/n  in  Fig.  3.  This 
results  in  uniaxial  beat  flow  toward  the  water- 
quenched  end,  and  reproducible  cooling  rates  that 
decrease  with  increasing  distance  away  from  the 
rapidly-cooled  quenched  end.  The  fast  cooling  rate 
at  the  quenched  end  causes  the  formation  of  hard, 
strong  martensite;  the  slower  cooled  end  transforms 
to  softer,  weaker  pearl ite  or  a  pearl ite  and  fer¬ 
rite  mixture,  depending  on  composition.  At  inter¬ 
mediate  locations  mixed  wrlcnxitc/poarlitc  micro- 
structures  result  frua  the  intermediate  cooling 
rates.  The  position  of  the  transition  from  marten¬ 
site  to  pearl itc  is  a  measure  of  the  hardening  res¬ 
ponse.  The  usual  way  of  assessing  the  microstruc¬ 
tural  gradient  along  the  length  of  an  cnd-quench 
bar  is  to  survey  the  Rockwell  C  hardness  measured 
on  flats  ground  along  the  side  of  the  end  quench 
test  bar.  This  is  possible  because  there  is  a 
pronounced  hardness  gradient  caused  by  the  micro- 
structural  gradient.  Although  surveying  hardness 
is  a  great  deal  easier  and  less  time  consuming 
than  directly  observing  the  microstructure  under 
a  microscope,  the  hardness  surveys  are  tedious 
because  up  to  60  individual  manual  hardness  meas¬ 
urements  may  be  required  to  survey  a  single  test 
bar,  a  procedure  taking  about  one  hour. 

Figure  4  compares  a  conventional  hardness 
survey  with  a  longitudinal  velocity  scan  of  an 
cnd-quench  test  bar  of  MSI  type  4140,  a  comnon, 
low  alloy,  heat- treatable  steel.  The  longitudinal 
acoustic  wave  velocity  was  measured  by  the  method 
described  above.)  The  acoustic  path  was  trans¬ 
verse  to  the  axis  of  the  test  bar.  Parallel  flats 
wvre  ground  on  opposite  sides  the  full  length  of 
the  test  bar;  the  acoustic  path  was  thus  along  bar 
diameters  through  the  thickness  of  the  bar  between 
the  ground  flats.  The  relative  acoustic  velocity 
is  plotted  as  Av/v^  .  the  fractional  change  in 
velocity  at  any  point  relative  to  Vq  the  velocity 


at  the  slower  cooled,  pcarlltic  end  of  the  test 
bar.  As  seen  in  the  plot,  the  velocity  decreased 
by  about  0.7  pet.  at  the  martensitic,  quenched  end. 
The  accuracy  of  the  relative  acoustic  velocity 
measurement  is  limited  by  the  uniformity  of  the 
thickness  and  the  deviation  from  perfect  parallel¬ 
ism  of  the  ground  flats.  This  is  estimated  by  be 
i  or  7  parts  in  10f'  which  is  about  21  of  the  actual 
range  of  velocities  measured.  Thus  the  sensitivity 
and  discrimination  of  the  velocity  measurement  is 
equal  to  or  better  than  that  of  the  hardness  meas¬ 
urements. 

DIAGRAM  or  A  JOMINV  TCST  IN  ACTION 


Fig.  3.  End-quench  test. 


Fig.  4.  End-quench  test  results  measured  in  type 
4140  steel.  The  two  curves  compare  the 
Microstructure  gradient  as  surveyed 
acoustically  and  by  conventional  hardness 
measurements. 

Figure  5  is  a  cross  plot  showing  the  relative 
velocity  change  Av/vn  ,  versus  the  hardness.  The 
correlation  is  virtually  linear  except  for  the 
extreme  values.  However,  it  is  the  intermediate 
values  that  are  most  important,  because  these 
occur  at  the  region  of  transition  from  martensite 
to  pearlite,  the  location  critical  to  the  measure? 
went  of  hardening  response. 
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Fig.  5.  Relative  acoustic  velcrity  vs.  hardness 
Rockwell  C  correlated  for  the  data  shown 
in  Fig.  A. 

The  computer  controlled  auto~atic  acoustic 
velocity  scan  requires  only  4  minutes  from  start  to 
finish  with  the  data  automatically  plotted  as 
Av/vq  vs.  position  in  the  bar.  This  contrasts 
favorably  with  the  hour  or  so  required  for  a  manual 
hardness  survey.  To  our  knowledge,  this  is  the 
first  time  an  end-quench  test  has  been  performed 
using  a  velocity  scan  rather  than  a  hardness  scan. 

In  addition  to  the  4140  steel  end-quench  test 
we  have  also  performed  similar  acoustic  velocity 
scons  on  end-quench  test  bars  of  types  52100,  4615, 
8640,  and  1096  steels  with  similarly  encouraging 
results.  V.’e  find  these  results  persuasive  that  • 
acoustic  velocity  measurements  can  be  used  to  map 
out  microstructural  variations  within  steel  objects. 

3.  Acoustic  Attenuation 

In  principle,  attenuation  measurements  are 
more  attractive  than  velocity  measurements  as  the 
basis  for  acoustic  HUE  microstructural  characteri¬ 
zation.  Where  velocity  is  weakly  a  function  of 
nicrostructure,  acoustic  attenuation  is  strongly 
affected  by  nicrostructure  because  of  scattering  at 
grain  boundaries,  second  phase  particles,  and  other 
microstructural  features.^  Moreover,  the  strong 
frequency  dependence  of  the  acoustic  attenuation 
coefficient  can  provide  additional  information 
related  to  nicrostructure. 

A  computer  interfaced  system  has  been  devel¬ 
oped  for  quickly  and  easily  measuring  the  attenu¬ 
ation  of  acoustic  waves  traveling  through  metal 
samples.  The  same  samples  used  for  the  velocity 
measurements  can  be  used  for  attenuation  measure¬ 
ment.  A  commercial  ultrasonic  transducer  launches 
longitudinal  waves  through  a  water  bath  at  normal 
incidence  to  the  sample  and  then  receives  the 
echoes. 

The  three- echo  method  is  used  to  determine 
the  attenuation,  as  described  by  Papadakis.^  This 
method  uses  the  measured  amplitudes  of  three  echoes 
to  solve  the  equations  for  these  echoes  for  any  of 
the  three  unknowns  in  the  system.  In  our  experi¬ 


ment,  these  unknowns  arc  the  pulse-echo  frequency 
response  of  the  system,  the  reflection  coefficient 
at  the  sample-water  interface,  and  the  attenuation 
through  the  sample.  Broadband  pulses  are  used  to 
obtain  information  over  the  frequency  range  of  an 
Octave  or  more.  The  front  face  eclio  and  first  two 
back  face  echoes  are  digitized  and  tiien  stored  in 
a  computer.  Since  the  dynamic  range  of  these 
echoes  can  be  40  dB  or  more,  many  samples  of  each 
echo  arc  averaged  to  improve  the  signal -to-noise 
ratio.  The  echoes  are  gated  and  separated  in  the 
computer,  and  the  moduli  of  their  Tourier  trans¬ 
forms  are  computed  by  an  FFT  program.  The  spectral 
components  of  thro  three  echoes  at  each  frequency 
arc  compared  in  the  three-echo  method  to  calculate 
the  attenuation  through  ttie  sample,  correcting  for 
transducer  response  and  also  for  reflection  and 
diffraction  losses.  Figure  6  illustrates  the 
Importance  of  including  the  diffraction  correction. 
The  scales  are  logarithmic,  log  dB/cm  versus  log  f, 
so  that  the  slopes  of  the  linear  plots  give  the 
powers  of  the  frequency  dependence.  These  curves 
are  plotted  by  a  computer  by  simply  connecting  128 
attenuation  data  points  within  the  passband.  No 
curve  fitting  techniques  have  been  used.  One 
measurement  produces  a  large  amount  of  attenuation 
information  about  the  sample. 


Fig.  6.  Attenuation  vs.  frequency  for  a  steel 
sample  illustrating  the  importance  of 
corrections  for  diffraction. 

Figure  7  shows  the  attenuation  in  two  differ¬ 
ent  steel  samples  with  identical  compositions  but 
different  microsctructurc.  The  curves  for  pearl ite 
and  martensite  show  frequency  dependences  close  to 
f**  indicating  that  Rayleigh  scattering  is  the 
dominant  mechanism  for  attenuation.  Pcarlitic 
data  was  obtained  from  two  different  transducers, 
as  shown  by  the  overlap  on  the  figure.  These 
results  are  comparable  to  previous  results  by 
Papadakis." 
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Fig.  7.  Attenuation  vs.  frequency  for  tv/o  specimens 
of  the  same  steel  heat  treated  differently 
to  have  martensite  in  one  but  pearl itc  in 
the  other. 

4.  Conclusions 

4.1  Absolute  acoustic  velocity  measurements  are 
subject  to  random  variations  that  limit  their 
utility  for  microstructural  NDE. 

4.2  Relative  velocity  change  measurements  can  be 
used  to  map  microstruclure  with  precision  in  a 
give'  piece  of  steel. 

4.3  Acoustic  attenuation  measurements  are  very 
sensitive  to  microstruclure,  and  they  have  good 
pote'tial  for  practical,  microstructural  NDE. 
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ABSTRACT 

The  measurement  of  acoustic  properties  can  be  used  for  the  nondestructive  characterisation  of  the 
mlcrostructure  of  materials.  We  have  measured  the  changes  In  longitudinal  acoustic  wave  velocity  and 
acoustic  attenuation  in  steel  specimens  whose  microstructure  and  properties  differ  widely  because  of 
differing  compositions  and  heat  treatment.  The  spatial  variation  ot  the  relative  acoustic  velocity  in 
standard  Jominy  cnd-qucnch  hardenahil tty  test  specimens  was  found  to  correlate  very  well  with  Rockwell  C 
hardness  scans,  indicating  a  potentially  practical  method  for  measuring  the  hardening  response  ot  heat 
treated  steel.  Absolute  velocity  measurements  on  steel  specimens  were  found  to  be  subject  to  random 
scatter  related  to  minor  compositional  variations;  this  limits  the  utility  of  absolute  velocity  measure 
moots  for  microstructural  NOT.  Attenuation  measurements  have  also  been  performed  on  steel  samples  with 
different  microstructures.  The  measurement  utilized  broadband  acoustic  pulses  corrected  for  transducer 
response,  liquid  buffer/solid  specimen  reflection,  and  diffraction  effects.  Attenuation  coefficients  were 
seen  to  be  proportional  to  frequency  squared  tor  martensite  and  to  the  fourth  power  ot  frequency  for 
pearlite.  Higher  attenuation  was  observed  for  pearlitic  than  for  martensitic  microstructures. 


Acoustic  methods  can  be  applied  to  the  nondes¬ 
tructive  characterization  of  the  microstructure  of 
materials.  This  paper  reports  the  initial  efforts 
of  our  Interdisciplinary  program  to  explore  micro- 
structural  acoustic  N01 .  The  experimental  materials 
were  steels. 

The  properties  of  materials,  especially  mechan¬ 
ical  properties  such  as  tensile  strength,  hardness, 
and  impact  toughness,  are  strongly  dependent  on 
microstructural  features  like  grain  size  and  shape, 
the  proportions  and  spatial  distribution  ot  the 
phases  present,  and  macroscale  heterogeneities  in 
mlcrostructure  resulting  from  local  compositional  or 
thermal-mechanical  processing  variations,  lliis  is 
true  of  both  metallic  and  ceramic  materials.  The 
microstructural  character  ot  materials  and  their 
related  properties  are  usually  assessed  by  micro¬ 
scopic  examination  or  by  direct  measurement  of  prop¬ 
erties;  such  procedures  are  often  costly,  and  are 
necessarily  destructive,  requiring  sacrificial  ex¬ 
ample  specimens. 

Steels  are  basically  alloys  of  Iron  and  carbon 
modified  by  minor  additions  of  other  elements. 

Steels  exhibit  remarkable  variability  of  mechanical 
properties.  These  useful  properties  are  dependent 
upon  the  rather  complex  mlcrostructure  of  steel 
which  can  be  closely  controlled  b.v  proper  selection 
of  composition  and  heat  treatment. 

The  dependence  of  steel  mlcrostructure  on  com- 
osltlon  and  heat  treatment  is  Illustrated  in  lig.  1. 
tgures  1(a)  and  1(b)  compare  the  microscopic  ap¬ 
pearance  of  polished  and  etched  steel  specimens  con¬ 
taining  O.i’  and  O.b  wl.  pet.  carbon  respectively; 
after  annealing  at  900*C  and  slow  cooling  in  air  to 
room  temperature  both  steels  consist  of  ferrite, 
the  white  grains,  and  pear  I  He,  the  darker  constit¬ 
uent.  lerrite  is  essentially  pure  iron;  pearlite 
Is  a  mixture  of  two  different  crystalline  phases 
(too  finely  dispersed  to  be  resolved  at  the  magnifi¬ 
cation  shown),  iron  carbide  (lofC)  and  lerrite. 

Since  virtually  all  the  carbon  resides  in  the  pearl¬ 
ite,  the  proportion  ot  pearlite  is  greater  In  the 
higher  carbon  steel,  figure  1(c)  shows  a  steel  with 
0.4  wt  pet  carbon  with  a  martensic  mlcrostructure 


rather  than  pearlite  and  ferrite.  This  steel  has 
been  cooled  rapidly  in  water  alter  annealing  at 
900't;  t  he'  fast  cooling  causes  the'  formation  of 
martensite,  a  supersaturated  solid  solution  ot  car¬ 
bon  in  iron.  Mechanical  properties  are  strongly 
affected  by  the  proportions  of  ferrite  and  pearlite 
and  by  the  existence  of  martensite  rather  than 
pearlite  and  ferrite.  Changes  in  mlcrostructure 
such  as  these,  enable  steels  to  be  heat  treated 
to, obtain  optimum  combinations  of  properties. 

AfOUSl  If  VliOCHY  Ml  AMIRI  MINI 

Our  measurements  cf  longitudinal  acoustic  wave 
velocity  In  metallic  samples  wore  made  with  a  com¬ 
puter-control!  system  developed  for  measuring 
acoustic  vtloc  fields  in  solid  samples  iinnersed 

in  a  liquid  hi  r.  The  measurements  use  a  two 

pulsp-ccho  ti'i  mque  that  cancels  out  the  effects 
ol  the  liquid  buffer,  lhr  measuring  method  is  il¬ 
lustrated  in  I ig.  >'.  Iwo  sets  of  tone  bursts 
several  rf  cycles  long  are  transmitted,  and  their 
reflections  from  the  specimen  are  received  by  a 
mechanically  scanning  commercial  transducer.  1  Ik' 
delay  between  the  two  pulses  is  adjusted  to  overlap 
the  back-face  echo  ol  the  first  pulse  with  the 
front-face  echo  cl  the  second  pulse,  and  this  sum  is 
gated  out  and  detected,  llio  product  term  is  then 
used  as  an  error  signal  in  a  phase-lock  loop,  which 
adjusts  the  frequency  to  keep  the  phase  difference 
between  the  two  echoes  constant.  I  ifectively,  the 
measurement  ol  phase  change  Introduced  by  propaga¬ 
tion  through  the  specimen  is  converted  to  that  ol 
frequency,  which  can  be  pcrlonned  with  great  pre¬ 
cision.  All  of  the  above  operations,  including  (lie 
mechanical  scanning  of  the  transducer,  data  collec¬ 
tion,  reduction,  and  display  are  controlled  by  a 
POP  11-.I4  minicomputer,  llio  resulting  system  pre¬ 
cision  is  about  1  part  in  10*'  ol  the  measured  longi¬ 
tudinal  acoustic  velocity,  lhr  accuracy  ol  the 
absolute  velocity  measurement  is  ,'  parts  in  10'*,  as 
dictated  by  the  accuracy  in  measuring  the  acoustic 
path  length,  the  specimen  thickness. 


(c) 


FIG.  1.  Different  steel  microstructures  resulting 
from  difference  in  composition  and  heat 
treatment.  1(a):  pearlite  and  ferrite  in 
normalized  0.2  wt  pet  C  steel; 
l(l>):  pearlite  and  ferrite  in  normalized 
0.6  wt  pet  C  steel;  1(c):  martensite  in 
water  quenched  0.4  wt  pet  C  steel. 
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FIG.  2.  Schematic  of  the  precision  acoustic 
velocity  measuring  system. 

ABSOl 01 1  VriOCilY  Mi  ASUWMi NTS 

Specimens  of  5  plain  carbon  steels  were  pre¬ 
pared  by  machining  1  cm  thick  flat  plates  oriented 
both  longitudinally  and  transversely  relative  to 
the  steel  bar  stock.  Since  the  acoustic  wave  prop¬ 
agation  direction  was  through  the  1  cm  plate  thick¬ 
ness,  acoustic  velocity  could  be  measured  both 
parallel  to  and  transverse  to  the  rolling  direction 
for  each  steel  composition  to  detect  the  influence 
of  any  preferred  orientation  or  crystal  texturing. 
The  compositions  of  the  steels  are  shown  in  Table  1. 
All  specimens  were  given  Identical  heat  treatment; 
they  were  heated  to  900'C  and  air  cooled  (normal¬ 
ized).  This  resulted  in  pearl ite/ferrite  micro- 
structures,  the  pearlite  varying  from  12  to  100 
volume  pet  over  the  range  of  carbon  contents. 

After  machining,  the  flat  specimens  were  lapped  to 
assure  that  their  flat  surfaces  were  parallel  to 
within  2,!>  pm, 

TABU  I 

Composition  of  Plain  Carbon  Steels  used  for 
Acoustic  Velocity  Measurements 


cel  Type 
(A LSI ). 

C 

Mn 

Si 

P 

s 

1010 

0.10 

0.45 

0.03 

0.007 

0.032 

1020 

0.20 

0.52 

0.21 

0.009 

0.07!! 

103  b 

0.33 

0.72 

0.17 

o.om 

0.075 

lObO 

0. 55 

0.00 

0.19 

0.012 

0.073 

1095 

0.93 

0.50 

0.22 

0.009 

0.030 

The  resulting  velocity  measurements  exhibited 
random  scatter  that  totally  obscured  the'  variation 
in  velocity  anticipated  from  the  variation  in 
microstructure,  Fig.  .1.  One  can  easily  calculate 
an  expected  acoustic  velocity  from  the^known  den¬ 
sities  of  pure  iron  ferrite*  and  le;C,*'  and  known 
elastic  constant  data  lor  ferrite/legC  mixtures.' 
Increasing  carbon  should  cause  a  decrease  in  the 
longitudinal  acoustic  velocity  linearly  propor¬ 
tional  to  carbon  content;  velocity  should  decrease 


0.85  pet  per  wt.  pet  carbon  In  the  steel,  all  other 
factors  being  identical.  The  scatter  observed  in 
the  aieasured  velocity  was  not  caused  by  variations 
In  preferred  orientation;  had  there  been  any  signi¬ 
ficant  degree  of  preferred  orientation,  the  longi¬ 
tudinal  and  transverse  velocities  would  differ  much 
more  than  the  slight  variations,  which  can  be  at¬ 
tributed  to  sample  thickness  variations,  scon  in 
Fig.  3.  We  believe  that  the  random  variations  in 
such  elements  as  manganese,  silicon,  sulfur  and 
others  always  present  in  at  least  trace  amounts, 
cause  variations  in  density  and  elastic  modulus 
and  thus  introduce  random  variations  in  acoustic 
vel  ocity  equal  to  or  greater  than  the  systematic 
effect  of  carbon  in  changing  the  microstructure. 

Tor  example,  we  have  calculated  that  the  range  of 
manganese  contents  in  our  specimens  (see  Table  I) 
causes  density  variations  large  enough  to  obscure 
all  the  above  change  attributable  to  variations  in 
microstructure.  Therefore  we  conclude  that  abso¬ 
lute  velocity  measurements  are  not  practical  means 
for  characterizing  steel  microstructure,  because 
random  compositional  variations  will  confuse  any 
attempt  to  calibrate  absolute  acoustic  velocity 
with  microstructure. 


MEASURED  ACOUSTIC  VELOCITY 


CARBON  (PERCENT ) 


FIG.  3.  Measured  acoustic  velocity  in  plain  carbon 
steels  of  varying  carbon  content.  The 
random  scatter  caused  by  minor  fluctuations 
In  alloying  elements  obscures  the  system¬ 
atic  variation  in  velocity  related  to  the 
microstructure. 

RELATIVE  VELOCITY  MEASUREMENTS 


By  relative  velocity  measurements  we  mean  mea¬ 
surements  of  acoustic  velocity  at  different  loca¬ 
tions  within  a  given  object.  If  the  material  compo¬ 
sition  is  uniform  within  the  object  (as  is  usually 
the  case),  variations  in  acoustic  velocity  can  be 
unequivocally  related  to  microstructura)  variations 
from  one  location  to  another.  There  are  many  situa¬ 
tions  where  microstructurc  is  not  uniform  and  the 
pattern  of  microstructura 1  variation  is  important, 
c.g.,  in  an  induction  hardened  steel  shaft  heat 
treated  to  bo  martensitic  on  the  outer  surface  but 
pearlitic  in  its  interior.  Velocity  variation  within 
•  given  object,  l.o.,  relative  velocity,  could  be 


used  to  map  out  the  microstructural  pattern.  Jominy 
end-quench  test  bars  provide  another  example  of 
steel  objects  containing  microstructural  gradients. 

We  successfully  used  relative  velocity  measurements 
to  survey  quantitatively  the  microstructurc  of  end 
quench  test  bars  of  several  different  alloy  steel 
compositions. 

The  Jominy  Test  or  End-Quench  Test  is  the 
standard  metallurgical  quali’y  control  procedure  for 
measuring  the  hardening  response  of  heat  treatable 
steels.  The  end-quench  test  specimen  is  a  one-inch 
diameter  round  bar,  four  inches  long,  lhe  bar  is 
heated  to  about  850'C.  It  is  then  placed  in  a  fix¬ 
ture  and  cooled  by  a  jet  of  cold  water  impinged  upon 
one  end  as  shown  in  Tig.  4.  This  results  in  uniaxial 
heat  flow  toward  the  water-quenched  end,  and  repro¬ 
ducible  cooling  rates  that  decrease  with  increasing 
distance  away  from  the  rapidly-cooled  quenched  end. 
The  fast  cooling  rate  at  the  quenched  end  causes  the 
formation  of  hard,  strong  martensite;  the  slower 
cooled  end  transforms  to  softer,  weaker  pearl ite  or 
a  pearlite  and  ferrite  mixture,  depending  on  compo¬ 
sition.  At  intermediate  locations  mixed  martensite/ 
pearlite  microstructures  result  from  the  intermedi¬ 
ate  cooling  rates.  The  position  of  the  transition 
from  martensite  to  pearlite  is  a  measure  of  the 
hardening  response.  The  usual  way  of  assessing  the 
microstructural  gradient  along  the  length  of  an  end- 
quench  bar  is  to  survey  the  Rockwell  C  hardness 
measured  on  flats  ground  along  the  side  of  the  end 
quench  test  bar.  This  is  possible  because  there  is 
a  pronounced  hardness  gradient  caused  by  the  micro- 
structural  gradient.  Although  surveying  hardness  is 
a  great  deal  easier  and  less  time  consuming  than  di¬ 
rectly  observing  the  microstructu^e  under  a  micro¬ 
scope,  the  hardness  surveys  are  tedious  because  up 
to  f>0  individual  manual  hardness  measurements  may  be 
required  to  survey  a  single  test  bar,  a  procedure 
taking  about  one  hour, 

DIAGRAM  OF  A  JOMINY  TEST  IN  ACTION 


FIG.  A.  End-quench  test. 

Figure  5  compares  a  conventional  hardness  survey 
with  a  longitudinal  velocity  scan  of  an  end-quench 
test  bar  of  AISl  type  4140.  a  common,  low  alloy, 
heat- treatable  stool.  The  longitudinal  acoustic  wave 
velocity  was  measured  by  the  method  described  above 
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and  frustrated  In  Fig.  ?.  The  acoustic  path  was 
transverse  to  the  axis  of  the  test  bar.  Parallel 
flats  were  ground  on  opposite  sides  the  full  length 
of  the  test  bar;  the  acoustic  path  was  thus  along 
bar  diameters  through  the  thickness  of  the  bar 
between  the  ground  tlats.  Ihe  relative  acoustic 
velocity  is  plotted  as  AV/Vq  ,  the  fractional 
change  in  velocity  at  any  point  relative  to  V0 
the  velocity  at  the  slower  cooled,  pearlitic  end  of 
the  test  bar.  As  seen  in  the  plot  the  velocity 
decreased  by  about  0.7  pet  at  the  martensitic, 
quenched  end.  Ihe  accuracy  of  the  relative  acous¬ 
tic  velocity  measurement  is  limited  by  the  uni¬ 
formity  of  the  thickness  and  the  deviation  from  per¬ 
fect  parallelism  of  the  ground  flats.  This  is  esti¬ 
mated  to  be  1  or  2  parts  in  lO^  which  is  about  2% 
of  the  actual  range  of  velocities  measured.  Thus 
the  sensitivity  and  discrimination  of  the  velocity 
measurement  is  equal  to  or  better  than  that  of  the 
hardness  measurements. 

Figure  6  Is  a  cross  plot  showing  the  relative 
velocity  change  AV/Vq  ,  versus  the  hardness.  The 
correlation  is  virtually  linear  except  for  the  ex¬ 
treme  values.  However,  it  is  the  intermediate 
values  that  are  most  important,  because  these  occur 
at  the  region  of  transition  from  martensite  to 
pearlitc,  the  location  critical  to  the  measurement 
of  hardening  response. 

The  computer  controlled  automatic  acoustic 
velocity  scan  requires  only  <1  minutes  from  start  to 
finish  with  the  data  automatically  plotted  as 
AV/V  vs  position  in  the  bar.  This  contrasts 
favorably  with  the  hour  or  so  required  for  a  manual 
hardness  survey.  To  our  knowledge,  this  is  the 
first  time  an  rnd-quench  test  has  been  performed 
using  a  velocity  scan  rather  than  a  hardness  scan. 

In  addition  to  the  4140  steel  cnd-quench  test 
we  have  also  performed  similar  acoustic  velocity 
scans  on  cnd-quench  test  bars  of  types  52100,  -1615, 
8640,  and  1095  steels  with  similarly  encouraging 
results.  We  find  these  results  persuasive  that 
acoustic  velocity  measurements  can  be  used  to  map 
out  microstructural  variations  within  steel  objects. 
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FJG.  5.  End-qucnch  test  results  measured  in  typo 
4140  steel.  Ihe  two  curves  compare  the 
micros  true  lure  gradient  as  surveyed 
acoustically  and  by  conventional  hardness 
measurements. 


FIG.  6.  Relative  acoustic  velocity  vs  hardness 

Rockwell  C  correlated  for  the  data  shown 
In  Fig.  5. 

ACOUSTIC  ATTINUAT 1  ON 

In  principle,  attenuation  measurements  are 
more  attractive  than  velocity  measurements  as  the 
basis  for  acoustic  NOE  microstructural  characteriza¬ 
tion.  Where  velocity  is  weakly  a  function  of  micro¬ 
structure,  acoustic  attenuation  is  strongly  affected 
by  microstructure  because  of  scattering  at  grain 
boundaries,  second  phase  particles  and  other  micro- 
structural  features. **  Moreover,  the  strong  fre¬ 
quency  dependence  of  the  acoustic  attenuation  coef¬ 
ficient  can  provide  additional  information  related 
to  microstructure. 

This  paper  reports  only  our  earliest,  prelimi¬ 
nary  efforts  to  characterize  microstructure  by  means 
of  acoustic  attenuation  measurements. 

A  computer  interfaced  system  has  been  developed 
for  quickly  and  easily  measuring  the  attenuation  ot 
acoustic  waves  traveling  through  metal  samples. 

The  same  samples  used  for  the  velocity  measurements 
can  be  used  for  attenuation  measurement.  A  commer¬ 
cial  ultrasonic  transducer  launches  longitudinal 
waves  through  a  water  hath  at  normal  incidence  to 
the  sample  and  then  receives  the  echoes.  Normal 
Incidence  is  attained  by  swiveling  the  transducer 
to  maximize  the  amplitude  of  the  front  lace  echo. 

A  narrowband  pulse  (i.e.,  a  tone  burst  with  at 
least  10  rf  cycles)  must  be  used  for  this  alignment 
because  the  various  components  of  a  broadband  pulse 
are  not  equally  affected  by  non-normal  incidence 
and  therefore  the  maximum  is  difficult  to  determine. 
Computer  controlled  digital  stepping  motors  can  also 
be  used  to  move  the  transducer  for  spatial  scanning. 

The  three-echo  method  is  used  to  determine  the 
attenuation,  as  described  by  Papadakis.*  Ibis 
method  uses  the  measured  amplitudes  of  three  echoes 
to  solve  the  equations  for  these  echoes  for  any  of 
the  three  unknowns  in  the  system.  In  our  experiment, 
those  unknowns  are  the  pulse-echo  frequency  response 
of  the  system,  the  reflection  coefficient  at  the 
sample-water  interface,  and  the  attenuation  through 
the  sample.  Broadband  pulses  arc  used  to  obtain 
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Information  over  the  frequency  ranqe  of  an  octave  or 
more,  The  front  face  echo  and  first  two  back  face 
echoes  are  digitized  and  then  stored  in  a  computer. 
Since  the  dynamic  range  of  these  echoes  can  be 
40  dll  or  more,  many  samples  of  each  echo  arc  aver¬ 
aged  to  improve  the  signal-to-noise  ratio.  The 
echoes  are  gated  and  separated  in  the  computer,  and 
the  moduli  of  their  Fourier  transforms  are  computed 
by  an  F F T  program.  Figure  7  shows  three  gated, 
averaged  puTses  with  their  frequency  spectra  super¬ 
imposed  for  a  typical  sample.  The  spectral  compo¬ 
nents  of  the  three  echoes  at  each  frequency  arc 
compared  in  the  three-echo  method  to  calculate  the 
attenuation  through  the  sample,  correcting  for 
transducer  response  and  also  for  reflection  and 
diffraction  losses.  Figure  8  illustrates  the  im¬ 
portance  of  including  the  diffraction  correction 
for  the  echoes  shown  in  Fig.  7.  The  scales  are 
logarithmic,  log  dB/cm  versus  log  f  ,  so  that  the 
slopes  of  the  linear  plots  give  the  powers  of  the 
frequency  dependence,  lhese  curves  are  plotted  by 
a  computer  and  an  analog  plotter  by  simply  connect¬ 
ing  128  attenuation  data  points  within  the  passband. 
No  curve  fitting  techniques  have  been  used.  One 
measurement  produces  a  large  amount  of  attenuation 
information  about  the  sample. 
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FIG.  8.  Attenuation  vs  frequency  for  a  steel 
sample  illustrating  the  importance  of 
corrections  for  diffraction. 

Figure  9  shows  the  attenuation  in  two  differ¬ 
ent  steel  samples  with  identical  composition  but 
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FIG.  7.  Three  echoes  of  the  same  broadband  pulse 
(plotted  in  both  the  time  and  frequency 
domains)  used  to  obtain  the  attenuation 
coefficient  in  the  specimen  for  the  range 
of  frequencies  in  the  pulse. 
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FIG.  9.  Attenuation  vs  frequency  for  two  specimens 
of  the  same  steel  heat  treated  differently 
to  have  martensite  In  one  but  pearlite  in 
the  other. 
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different  microstructure.  lhe  curve  for  pearlite 
shows  a  frequency  dependence  very  close  to  fA 
Indicating  that  Rayleigh  scattering  Is  the  dominant 
Mechanism  for  attenuation,  lhe  curve  for  martensite 
exhibits  lower  values  oi  attenuation,  with  a  fre¬ 
quency  dependence  nearer  to  f ,  Indicating  a  dif¬ 
ferent  mechanism  for  attenuation  from  that  In  mar¬ 
tensite.  lhese  results  are  comparable  to  previous 
results  by  Papadakis.5* 

CONCLUSIONS 

1.  Absolute  acoustic  velocity  measurements 
are  subject  to  random  variations  that  limit  their 
utility  for  microstructural  NHL. 

?.  Relative  velocity  change  measurements  can 
be  used  to  map  microstructure  with  precision  in  a 
given  piece  of  steel. 

3.  Acoustic  attenuation  measurements  are  very 
sensitive  to  microstructure,  and  they  have  good 
potential  for  practical,  microstructural  NOE. 
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Abstract 

An  automatic  system  is  described  for  measuring 
two-dimensional  acoustic  velocity  fields  in  solid 
samples.  The  measurement  is  performed  by  a  com¬ 
puter-controlled,  mechanically-scanned  transducer 
in  a  liquid  bath  and  is  based  on  measuring  the 
phase  delay  of  the  acoustic  wave  by  a  two-pulse 
echo  method.  Applications  include  measuring  stress 
fields  due  both  to  externally-applied  and  residual 
stresses,  and  raicrostructure  studies  of  solid 
samples. 


1.  Introduction 

Recently,  there  has  been  a  resurgence  of 
interest  in  studies  of  material  properties  by 
acoustic  waves.  Two  important  reasons  for  that 
development  include  the  availability  of  dedicated 
computers  for  control  of  experiments  and  massive 
data  processing,  and  the  important  expansion  of  the 
field  of  nondestructive  materials  evaluation  (NOE), 
We  report  here  on  the  development  and  the  applica¬ 
tion  to  h'DE  of  a  measuring  system  which  automati¬ 
cally  produces  two-dimensional  images  of  stress 
distributions  in  metals  and  other  solid  materials. 

The  measuring  technique  is  based  on  the  two- 
pulsc  echo  method  developed  some  time  ago,  and  is 
applied  normally  to  measuring  propagation  velocity 
in  solid  samples  with  bonded  ultrasonic  trans¬ 
ducers.2  By  contrast,  our  experiment  is  carried 
out  with  a  liquid  buffer  so  that  samples  of  a 
variety  of  shapes,  sizes,  and  compositions  can  be 
tested  |>y  a  mechanically- scanned  transducer.  The 
entire  two-dimensional  scanning  operation  is  com¬ 
pletely  automated  and  computer-controlled,  so  that 
complete  scans  containing  several  hundred  data 
points  can  be  made  and  the  results  displayed  in 
about  thirty  minutes.  Vie  have  already  reported  on 
an  early,  manually-operated  version  of  the  instru¬ 
ment,  which  did  not  have  the  feature  of  cancelling 
out  the  water  buffer  propagation  effects,  and  . 
provided  only  preliminary  two-dimensional  scans. 

A  norc  recent  publication  has  reported  some  results 
of  acousto-elastic  measurements  in  metals,  per¬ 
formed  on  the  present  system.  We  present  here  the 


details  of  the  measuring  system,  and  then  discuss 
some  additional  applications. 

2.  Measuring  System 

The  measuring  system  consists  of  a  liquid 
filled  tank  which  is  constructed  as  part  of  a 
mechanical  rig  for  applying  calibrated  tensile  or 
compressive  loads  to  the  plane-paral lei  specimen. 

As  the  acoustic  wave  velocity  is  changed  by  the 
applied  stress,  the  variation  of  stress  over  the 
cross-section  of  a  sample  induced  by  external 
loading  can  be  measuring  by  a  mechanically-scanned 
transducer  while  the  load  is  being  applied.  The 
homemade,  hydra ul ical ly-operated  loading  rig  can 
apply  loads  up  to  -5  x  lCrN  . 

A  piston  transducer  is  scanned  mechanically 
in  the  vertical  and  horizontal  directions  in  the 
liquid-filled  tank  by  a  pair  of  computer-control  led 
stepping  motors.  The  liquid  employed  is  either 
water  or  ethylene  glycol,  the  latter  being  more 
resistant  to  metal  sample  corrosion,  the  build-uo 
of  mineral  deposits,  and  the  occurrence  of  bubbles 
at  the  sample  surface.  Typically,  an  unfocused 
transducer  is  used  to  launch  a  beam  of  longitudinal 
acoustic  waves  of  about  3  mm  diameter  at  around 
J2  MHz. 

The  measurement  principle  is  illustrated  in 
Fig.  J.  Two  rf  tone  bursts  a  few  us  long  are 
applied  to  the  transducer,  with  a  variable  time 
delay  between  them.  At  .  The  corresponding  acous¬ 
tic  waves  are  launched  such  that  they  impinge  on 
the  sample  at  normal  incidence,  and  longitudinal 
waves  only  are  excited  In  the  sample.  At  the 
liquid  buffer/solid  sample  interface,  the  waves  are 
partially  transmitted  and  reflected;  the  reflected 
part  propagates  back  and  is  received  by  the  trans¬ 
ducer,  while  the  transmitted  part  of  the  acoustic 
energy  undergoes  further  reflections  at  the  sample 
back  and  front  faces  as  it  bounces  back  and  forth 
inside  the  sample.  Each  one  of  these  transits 
Toads  to  some  acoustic  energy  being  transmitted 
back  to  the  transducer,  resulting  in  a  correspond¬ 
ing  received  set  of  rf  pulses.  Because  a  liquid 
path  is  employed,  the  reflection  coefficient,  both 
at  the  front  and  the  back  of  the  sample,  is  not 
dependent  on  the  bonding  material  and  should  be 
invariant  with  the  transducer  position. 
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Fig.  1.  the  principle  of  the  tvio-pulse  echo  measuring 
technique  in  a  liquid  buffer 

In  order  to  determine  the  local  velocity  of  longi¬ 
tudinal  wave  propagation  in  the  sample.  At  is  adjus¬ 
ted  so  that  the  received  rf  pulses,  due  to  the  front- 
face  echo  of  the  second  pulse  and  the  back-face  of  the 
first  pulse,  overlap.  This  sum  is  then  qated  out  from 
•11  the  other  received  signals,  and  the  phase  shift 
between  the  two  echoes  is  measured,  corresponding  to  a 
double  transit  of  the  wave  propagation  throuqh  the 
sample.  Thus,  by  usinq  the  front-face  echo  as  the 
reference  pulse,  the  wave  propagation  effects  through 
the  water  buffer  and  the  possible  additional  phase 
shifts  in  the  outside  circuitry0  are  eliminated  from 
this  measurement. 


3.  Electronic  Circuitry 
3.1  Relative  Velocity 

The  requirement  for  a  very  precise  measurement  of 
the  phase  difference  between  the  front-face  and  the 
back-face  echoes,  which  could  at  the  same  time  be  auto¬ 
mated,  was  met  by  an  indirect  measurement,  based  on  a 
nulling  technique.'  Neglecting  the  phase  shift  due  to 
the  reflection  at  the  sample  back-face,  the  phase 
change  which  the  acoustic  signal  undergoes  in  propagat¬ 
ing  through  the  sample  of  thickness  d  and  back  again 
to  the  front  face  is 


where  f  and  v  arc  the  wave  frequency  and  velocity, 
respectively.  For  reconstructing  velocity  profiles  in 
samples,  we  need  only  be  concerned  with  velocity  or 
transit  time  changes  in  moving  from  one  point  to  the 
next.  If  wc  vary  the  carrier  frequency  of  the  rf 
pulses  by  such  amounts  as  to  keep  $  -  const  at  every 
point  during  a  scan,  we  obtain  from  Eq.  (1) 


Af  Av .  Ad 
f  *  v  *  d  ' 


(2) 


It  has  been  shown  that  when  stress  is  applied, 
both  the  relative  change  in  velocity  Av/v  and  the 
relative  change  in  thickness  Ad/d  are  proportional 
to  the  applied  stress.**  Thus,  the  contours  of  con¬ 
stant  stress  in  a  sample  arc  eouivalent  to  contours  of 
constant  frequency  for  which  $  =  const  ,  It  is  also 
seen  that  the  constant  phase  shifts  due  to  reflections 
drop  out  of  this  relative  measurement,  hut  arc  needed 
for  an  absolute  measurement  of  velocity. 


lent  and  Is  highly  amplitude-sensitive.  We  therefore 
adopted  a  sampling  technique  in  which  the  basic  pulse 
repetition  rale  is  0  kHz  but  one  of  the  pulses  has  its 
sign  reversed  at  a  1  kilt  rate.  This  makes  it  possible 
to  use  a  square  law  detector  to  obtain  a  1  kilt  output 
which  depends  only  on  the  product  of  the  two  pulses  of 
Interest,  lhis  product  term 'in  turn  has  zero  ampli¬ 
tude  when  the  two  rf  pulses  have  carriers  n/2  out 
of  phase  with  each  other.  A  lock-in  amplifier  is 
employed  to  compare  the  1  kHz  detected  signal  with  a 
1  kHz  reference  so  that  the  system  can  be  made 
extremely  sensitive.  Furthermore,  the  output  of  the 
lock-in  amplifier  is  itself  employed  to  control  a 
phase-lock  loop,  which  changes  the  frequency  of  the 
Input  carrier  so  as  to  set  the  phase  difference  bet¬ 
ween  the  carriers  of  the  two  return  pulses  at  r/2  . 
Thus,  the  measurement  system  is  entirely  automatic. 

The  circuit  for  performing  the  measurement  and 
producing  results  in  the  form  of  hard-copy  contour 
plots  or  pictorial  scan  converter  output  is  shown  in 
Fig.  2.  A  free-running  oscillator  at  8  kHz  sets  the 
square  wave  pulse  repetition  rate  and  triggers  the 
excitation  of  the  sets  of  two-pulse  sequences.  The 
second  pulse  is  further  modulated  by  a  1  kHz  signal, 
applied  through  a  transistor  switch  (T S )  in  such  a 
way  that  at  the  8  kHz  repetition  rate,  the  square 
waves  change  sign  after  every  four  pulses,  the  two 
sets  of  square  waves  are  then  mired  with  the  cutout  of 
a  free-running,  highly  stable  synthesizer  (Hewlett- 
Packard  8660C);  the  resulting  rf  tone  bursts  are  added, 
amplified  by  a  40  dB  power  amplifier  and  transmitted 
by  the  ultrasonic  transducer. 


Fig.  2.  Detailed  block  diagram  of  the  measuring 

system.  PC,  -  pulse  generators;  VD  -  variable 
delay;  PA  -  power  amplifier;  RA  -  limiter 
and  receiving  amplifier;  SCOPE  -  oscilloscope 
([ton i tor;  LOCK-IN  -  lock  in  amplifier;  SYNTH  - 
synthesizer  of  rf  carrier. 


The  received  signal  from  the  same  transducer  is 
passed  through  a  limiter  which  protects  the  sensitive 
receiver  from  the  strong  transmitted  signal,  and  then 
amplified  by  a  CO  dR  variable  gain  receiving  amplifier. 
That  signal  is  displayed  by  the  oscilloscope,  at  which 
time  various  adjustments  and  calibrations  are  per¬ 
formed.  Only  the  desired  overlap  of  the  front  and 
back-face  echoes  indicated  in  fig.  1  is  then  gated 
out  and  detected  by  a  balanced  detector  (the  square 
law  detector).  The  siqnal  component  at  1  kHz  is  then 
synchronously  detected  with  a  lock-in  amplifier  with 
a  tuned  input  signal  channel. 


The  system  is  required  to  measure  the  phase  dif¬ 
ference  between  the  carriers  of  two  coincident  pulses. 
This  could  be  done  by  varying  the  amplitudes  of  the 
two  pulses  so  that  they  tend  to  cancel  when  they  are 
«  out  of  phase,  but  such  a  technique  Is  not  conven- 


ln  analyzing  the  received  and  gated  siqnals,  wo 
may  take  the  contribution  from  the  first  rf  pulse 
signal  to  be  of  the  form 

A(t)  cos  (2xft  ♦  ♦)  , 


(3) 


where  A(t)  Is  the  pulse  modulation  at  a  repetition  these  calculations  Instead  of  the  unknown  v  ,  a  value 
rate  of  8  kHz  .  Similarly,  the  contribution  from  the  very  close  to  but  not  exactly  equal  to  an  Integer  Is 

Second  pulse,  which  has  not  propagated  through  the  obtained  for  n  .  That  number  Is  then  rounded  off  to 

sample,  has  the  form  its  Integer  value,  and  v  Is  calculated  from  fq.  (8). 

An  Iterative  scheme  of  several  steps  could  also  be  set 

A(t)B(t)  cos  2wft  ,  (4)  up.  but  this  first  iteration  gives  us  good  accuracy. 


where  B(t)  is  the  1  kilt  square  wave  modulation  of 
the  second  pulse.  Vie  note  that  B(t)  •  il  and  that 
there  arc  eight  samples  (pulses)  per  cycle  of  B(t)  . 
If  we  now  square  the  sum  of  Eqs.  (3)  and  (4)  and  take 
only  the  contribution  at  1  kHz,  the  corresponding 
tnput  to  the  lock-in  amplifier  becomes 

A?B  cos  4  .  (5) 

The  output  is  then  a  DC  signal  proportional  to  R  , 
where 


3.3  Diffraction  Correction 

The  phase  change  term  due  to  diffraction  4j  in 
Eq.  (8),  was  derived  in  an  approximate  form,  valid  in 
the  paraxial  limit.  The  expression  is  very  suitable 
for  computer  calculation,  and  the  reader  is  referred 
for  the  derivation  Steps  leading  to  it  to  a  forth¬ 
coming  publication." 

We  are  Interested  in  the  detected  pressure  across 
the  face  of  a  receiver  of  radius  a  at  distance  z 
away.  The  required  expression  is 


R  *  A  cos  4  cos  4  ,  (6) 

where  4  is  a  constant  phase  shift  between  the  refer¬ 
ence  and  signal  paths  of  the  1  kHz  modulating  signal, 
and  4  is  the  unknown  desired  phase  shift. 

The  DC  signal  given  by  Eq.  (6)  is  then  used  as  an 
error  signal  in  a  phase-lock  loop  which  varies  the 
synthesizer  output  frequency  in  such  a  way  that  R  1  0 

or 
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Now  we  introduce  the  approximation,  k  »  o  , 


4  -  (2n  *  3)n/2  ,  (7) 

since  the  loop  locks  onto  every  second  null  of  cos  4  • 
The  lock-in  amplifier  output  is  first  integrated  (the 
lock-in  amplifier  integration  constant  is  1  ms,  that 
of  the  integrator  is  10  ms),  and  then  applied  to  the 
FM  input  of  the  synthesizer,  so  that  the  resulting 
frequency  f  adjusts  automatically  to  keep  R  -  0 
In  Eq.  (6).  This  frequency  is  then  read  into  the 
co-puter  via  an  electronic  counter,  and  the  single¬ 
point  measurement  is  completed. 

3.2  Absolute  Velocity 

For  performing  absolute  velocity  measurements, 
the  rf  carrier  frequency  is  varied  over  a  wide  range, 
and  several  neighboring  nulls  corresponding  to  Eq.  (7) 
are  recorded.  We  write 

4nf  d 

4t  -  — +  *  -  4d  -  (2n  +  3)ir/2  ,  (8) 

where  4t  is  the  total  phase  shift  of  the  back-face 
echo,  taking  into  account  the  n  shift  from  the 
sarple-1 iquid  interface,  and  4d  is  the  contribution 
due  to  diffraction." 


and  substituting  S  =  Xz/a^  ,  where  X  =  2n/k  ,  y  =  aa  , 
Eq.  (10)  becomes 


*7.  y 


This  is  the  final  expression  used  in  computing 
the  phase  correction  due  to  diffraction,  i.e.,  we  take 


This  expression  can  be  evaluated  very  easily  on 
the  computer,  and  that  is  what  we  do  in  performing  the 
absolute  velocity  measurements.  As  an  illustration, 
the  results  of  calculating  Eq.  (13)  are  plotted  in 
Fig.  3,  both  for  the  phase  and  the  attenuation  correc¬ 
tions  (the  latter  is  important  for  correcting  attenu¬ 
ation  measurements).  These  results  agree  well  with 
the  more  elaborate  calculations  reported  in  the  liter¬ 
ature.5*' 


4.  Measuring  Procedure  and  Computer  Interfacing 


An  interactive  computer  program  has  beep  written 
which  solves  for  v  in  two  steps.  First,  from  the 
slope  of  f n+k  vs.  k  line  for  several  measured 
fn4)t  ,  an  approximate  value  of  velocity  v  is  calcu¬ 
lated  from  the  formula 

*1  ■  -  *„>/*  .  (9) 

This  Is  not  accurate  enough  for  the  final  calculations 
because  it  requires  taking  the  difference  of  two  large 
quantities. 

The  value  Vj  is  then  used  for  generating  the 
diffraction  phase  correction,  calculated  numerically 
for  our  acoustic  beam  diameter,  samnlc  thickness,  and 
transducer-sample  distance.  Eq.  (8)  is  next  solved 
for  n  using  v  «  Vj  ;  since  Vj  is  used  in  all 


Initial  calibration  consists  of  mechanically 
aligning  the  transducer  with  respect  to  the  sample  for 
maximum  front-face  reflection,  and  then  adjusting 
separately  the  pulse  widths  of  Pf,l  and  PC, 2  (refer  to 
Fig.  2)  such  that  echoes  from  one  pulse  do  not  qverlap, 
depending  on  sample  thickness  and  velocity  of  propaga¬ 
tion.  Tor  typical  samples  of  1  cm  thick  aluminum  or 
steel,  the  required  pulse  widths  range  from  1  to  sev¬ 
eral  ps  .  Next,  the  delay  of  the  second  pulse  (VD) 
is  adjusted  for  total  overlap  of  the  front-face  echo 
of  the  second  pulse  with  the  back-face  echo  of  the 
first  pulse  (a  few  ps  ),  and  the  attenuation  (VA) 
adjusted  for  the  two  pulses  to  have  about  equal  ampli¬ 
tudes.  Finally,  the  delay  of  the  gating  signal 
and  the  appropriate  pulse  width  of  are  adjusted 
to  pick  out  the  desired  sum,  and  the  center  frequency 
and  the  FM  level  of  the  synthesizer  are  adjusted  for 
appropriate  excursions  for  proper  frequency  null 


49  - 


tracking.  Typically,  for  the  above  stated  conditions, 
the  I'M  level  Is  -50  kite,  where  Af  -  300  kite  corres¬ 
ponds  to  A'k  •  ?*  at  f  -  12  Mite  .  The  measurement  _g 
precision  routinely  corresponds  to  A v/v  of  about  10 
with  r.otal  samples  of  the  order  of  1  cm  thick.  This 
ts  about  the  same  or  better  than  that  obtained  in  the 
non- scanned,  bonded  transducer  experiments  based  on 
pulse  superposition,  or  sing-around  methods. ‘ 


Fig.  3.  Computed  diffraction  corrections  to  phase 
and  amplitude 

The  measuring  procedure  for  scanning  a  sample 
Involves  defining  the  grid  of  transducer  positions  in 
computer  software,  and  after  the  initial  sample  and 
transducer  alignment  and  electronics  calibration,  a 
corresponding  array  of  frequency  values  of  nulls  of 
Eq.  (6)  is  generated  under  complete  computer  control, 
at  the  rate  of  a  few  seconds  per  point.  The  data 
points  are  displayed  either  as  contour  plots  of  equal 
velocities  or  as  equivalent  scan  converter  output.*1 

5,  Applications 

The  system  we  have  described  is  routinely  used  by 
a  varied  group  of  electrical  and  mechanical  engineers 
and  materials  scientists.  We  have  already  reported 
some  exciting  measurements  of  third  order  elastic  con¬ 
stants  and  comparisons  between  measured  and  calculated 
stress  fields  in  externally-loaded  samples  of  differ¬ 
ent  shapes.  Such  contour  plots  can  also  give  us  very 
useful  qualitative  information  about  stresses  around 
cracks  and  can  be  used  to  evaluate  crack  stress  inten¬ 
sity  factors.  More  recent  work  has  included  similar 
measurements  of  residual  stress  fields  in  samples  with 
stresses  built  into  them  by  a  process  of  hydrostatic 
extrusion.  Absolute  values  of  radial  and  hoop  stresses 
can  then  be  obtained  by  suitable  integration  of  contour 
data.  We  arc  just  at  the  beginning  of  these  latter, 
quantitative  applications,  and  the  prospects  for  NDE 
are  exciting  indeed. 

Ke  present  here  a  recent  example  of  the  applica¬ 
tion  of  this  system  to  the  study  of  a  standard  end- 
quench  Joniny  bar,  which  is  a  steel  rod  with  a  widely 
varying  hardness  along  its  axis.  The  rod  is  normally 
used  as  a  standard  for  measuring  the  hardening  response 
to  heat  treatment  of  steel  samples.  It  is  prepared  by 
Introducing  a  temperature  gradient  along  the  bar  axis, 
so  that  the  resultant  gradient  In  cooling  rate  causes 
a  marked  gradation  in  steel  mtcroslructure  and  conse¬ 
quently  hardness.  Fig.  <1  shows  the  results  of  hard¬ 
ness  (Rockwell  C)  measurements  together  with  a  one¬ 


dimensional  acoustic  velocity  scan,  showing  strong  cor¬ 
relation  between  the  two  curves.  It  Is  seen  that  the 
velocity  of  propagation  decreases  as  hardness  is 
Increased. 


DISTANCE  FROM  OUENCMEO  END  OF 
SPECIMEN  1  mm  ) 


Fig.  4.  Spatial  dependences  of  Rockwell  C  hardness 

(points)  and  ultrasonic  velocity  (line)  of  an 
end-quench  Jominy  bar  of  A1S1  4140  steel. 

V  «=  5.9  1jii/s 

Absolute  velocity  measurements  are  also  being 
carried  out  with  this  system,  as  described  in  the  pre¬ 
vious  section.  This  method  is  not  as  accurate  as  the 
relative  measurement  described,  because  measurements 
over  a  wide  range  of  frequencies  must  be  made,  and 
extraneous  factors  such  as  diffraction  corrections, 
amounting  to  about  0.15"  in  our  case,  should  be  taken 
into  account. ^  We  are  also  limited  by  the  thickness 
accuracy  of  our  samples  to  about  2  um  on  1  cm  samples, 
which  limits  the  ultimate  total  accuracy  to  about 
2  x  10’4  .  This  is  consistent  with  the  general  state¬ 
ment  that  the  precision  of  relative  velocity  scans  is 
about  two  orders  of  magnitude  better  than  the  accuracy 
of  absolute  velocity  measurements.^  However,  if  thick¬ 
ness  could  be  measured  accurately,  the  precision  of  the 
absolute  velocity  measurement  would  be  improved  accor¬ 
dingly. 

Finally,  we  are  now  preparing  a  shear  wave  scanner 
with  specially  constructed  contact  transducers  to  be 
used  in  conjunction  with  this  system.  Shear  wave  data 
will  enable  us  to  get  the  values  of  both  principal 
stresses  individually,  by  rotating  the  transducer  with 
respect  to  the  axis  of  preferred  orientation,  instead 
of  just  the  sum  of  the  principal  stresses  produced  by 
the  above  longitudinal  wave  data.1* 
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